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Abstract 

This is a pedagogical review on primordial non-Gaussianities from inflation models. We 
introduce formalisms and techniques that are used to compute such quantities. We review 
different mechanisms which can generate observable large non-Gaussianities during inflation, 
and distinctive signatures they leave on the non-Gaussian profiles. They are potentially 
powerful probes to the dynamics of inflation. We also provide a non-technical and qualitative 
summary of the main results and underlying physics. 
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1 Introduction 



An ambitious goal of modern cosmology is to understand the origin of our Universe, all the 
way to its very beginning. To what extend this can be achieved largely depends on what 
type of observational data we are able to get. Thanks to many modern experiments, we are 
really making progress in this direction. 

One of the representative experiments is the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite [THS]. It measures the light coming from the last scattering surface about 
13.7 billions years ago. This cosmic microwave background (CMB) is emitted at about 
379,000 years after the Big Bang, when electrons and protons combine to form neutral 
hydrogen atoms and photons start to travel freely through the space. Our Universe was 
very young at that moment and the large scale fluctuations were still developing at linear 
level. So the CMB actually carries valuable information much earlier than itself, which can 
potentially tell us about the origin of the Big Bang. 

There are two amazing facts about the CMB temperature map. On the one hand, it is 
extremely isotropic, despite the fact that the causally connected region at the time when 
CMB formed spans an angle of only about 0.8 degree on the sky today. On the other hand, 
we do observe small fluctuations, with AT/T ~ 10~ 5 . 

The inflationary scenario [?H9] naturally solves the above two puzzles. It was proposed 
nearly thirty years ago to address some of the basic problems of the Big Bang cosmology, 
namely, why the universe is so homogeneous and isotropic. In this scenario, our universe was 
once dominated by dark energy and had gone through an accelerated expansion phase, during 
which a Hubble size patch was stretched by more than 60 efolds or so. Inhomogeneities and 
large curvature were stretched away by this inflationary epoch, making our current observable 
universe very homogenous and flat. In the mean while, the fields that were responsible for 
and participated in this inflationary phase did have quantum fluctuations. These fluctuations 
also got stretched and imprinted at superhorizon scales. Later they reentered the horizon 
and seeded the large scale structures today [TOHbl] . 

The inflationary scenario has several generic predictions on the properties of the density 
perturbations that seed the large scale structures: 

• They are primordial. Namely, they were laid down at superhorizon scales and entering 
the horizon after the Big Bang. 

• They are approximately scale-invariant. This is because, during the 60 efolds, each 
mode experiences the similar expansion when they are stretched across the horizon. 

• They are approximately Gaussian. In simplest slow-roll inflation models, the inflaton 
is freely propagating in the inflationary background at the leading order. This is also 
found to be true in most of the other models and for different inflationary mechanisms. 
So the tiny primordial fluctuations can be treated as nearly Gaussian. 
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The CMB temperature anisotropy is the ideal data that we can use to test these pre- 
dictions. The obvious first step is to analyze their two-point correlation functions, i.e. the 
power spectrum. All the above predictions are verified to some extent [lj. The presence 
of the baryon acoustic oscillations proves that the density perturbations are indeed present 
at the superhorizon scales and reentering the horizon as the horizon expands after the Big 
Bang. The spectral index, n s = 0.963 ± 0.012, is very close to one, therefore, the den- 
sity perturbations are nearly scale invariant. Several generic types of non-Gaussianities are 
constrained to be less than one thousandth of the leading Gaussian component. 

But is this enough? 

Experimentally, the amplitude and the scale-dependence of the power spectrum consist 
of about 1000 numbers for WMAP. For the Planck satellite, this will be increased up to 
about 3000. However, we have about one million pixels in the WMAP temperature map 
alone, and 60 millions for Planck. So the information that we got so far is highly compressed 
comparing to what the data could offer in principle. This high compression is only justified if 
the density perturbations are Gaussian within the ultimate sensitivities of our experiments, 
so all the properties of the map is determined by the two-point function. Otherwise, we are 
expecting a lot more information in the non-Gaussian components. 

Theoretically, inflation still remains as a paradigm. We do not know what kind of fields 
are responsible for the inflation. We do not know their Lagrangian. We also would like to 
distinguish inflation from other alternatives. Being our very first data on quantum gravity, we 
would like to extract the maximum number of information from the CMB map to understand 
aspects of the quantum gravity. All these motivate us to go beyond the power spectrum. 

To give an analogy, in particle physics, two-point correlation functions of fields describe 
freely propagating particles in Minkowski spacetime. More interesting objects are their 
higher order correlations. Measuring these are the goals of particle colliders. Similarly, 
the power spectrum here describes the freely propagating particles in the inflationary back- 
ground. To find out more about their interaction details and break the degeneracies among 
models, we need higher order correlation functions, namely non-Gaussianities. So the role 
non-Gaussianities play for the very early universe is similar to the role colliders play for 
particle physics. 

With these motivations in mind, in this review, we explore various mechanisms that 
can generate potentially observable primordial non-Gaussianities, and are consistent with 
the current results of power spectrum. We will not take the approach of reviewing models 
one by one. Rather, we divide them into different categories, such that models in each 
category share the same physical aspect which leaves a unique fingerprint on primordial 
non-Gaussianities. On the one hand, if any such non-Gaussianity is observed, we know what 
we have learned concretely in terms of fundamental physics; on the other hand, explicit 
forms of non-Gaussianities resulted from this exploration provide important clues on how 
they should be searched in data. Even if the primordial density perturbations were perfectly 
Gaussian, to test it, we would still go through these analyses until various well- motivated 
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non-Gaussian forms are properly constrained. 
1.1 Road map 

The following is the outline of the review. For readers who would like to get a quick and 
qualitative understanding of the main results instead of technical details, we also provide a 
shortcut after the outline. 

In Sec. El we review the essential features of the inflation model and density perturbations. 

In Sec. El we review the first-principle in-in formalism and related techniques that will 
be used to calculate the correlation functions in time-dependent background. 

In Sec. HI we compute the scalar three-point function in the simplest slow-roll model. We 
list the essential assumptions that lead to the conclusion that the non-Gaussianities in this 
model is too small to be observed. 

In Sec. 0, we review aspects of inflation model building, emphasizing various generic 
problems which suggest that the realistic model may not be the algebraically simplest. We 
also introduce some terminologies used to describe properties of non-Gaussianities. 

Sec. El [7J and M contain the main results of this review. We study various mechanisms 
that can lead to large non-Gaussianities, and their distinctive predictions in terms of the 
non- Gaussian profile. 

In Sec. [91 we give a qualitative summary of the main results in this review. Before 
conclusion, we discuss several useful relations among different non-Gaussianities. 

Here is a road map for readers who wish to have a non-technical explanation and un- 
derstanding of our main results. After reading the short review on the inflation model and 
density perturbations in Sec. [H one may read the first and the last paragraph of Sec. @] to 
get an idea of the no-go statement, and then directly proceed to read Sec. [5j After that, one 
may jump to Sec. [9] where the main results are summarized in non-technical terms. 

The subject of the primordial non-Gaussianities is a fast-growing one. There exists many 
nice reviews and books in this and closely related subjects. The introductions to inflation 
and density perturbations can be found in many textbooks [T5T - KU] and reviews [2~TH2"o] . In- 
flationary model buildings in particle physics, supergravity and string theory are reviewed 
in Ref. [26H52] . Comprehensive reviews on the developments of theories and observations of 
primordial non-Gaussianities up to mid 2004 can be found in Ref. [55fl54"] . Recent comprehen- 
sive reviews on theoretical and observational developments on the bispectrum detection in 
CMB and large scale structure has appeared in Ref. [55115(3] . A recent comprehensive review 
on non-Gaussianities from the second order post inflationary evolution of CMB, which acts 
as contaminations of the primordial non-Gaussianities, has appeared in Ref. [57]. A recent 
review on how primordial non-Gaussianities can be generated in alternatives to inflation can 
be found in Ref. [58] . 
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2 Inflation and density perturbations 



In this section, we give a quick review on basic elements of inflation and density perturba- 
tions. We consider the simplest slow-roll inflation. The action is 



S 



d A Xyf—gR + / d 4 Xy/—g 



1 



<Td M <R0 - V((f>) 



(2.1) 



The first term is the Einstein-Hilbert action. The second term describes a canonical scalar 
field coupled to gravity through the metric g^ v . This is the inflaton, which stays on the 
potential V(<p) and creates the vacuum energy that drives the inflation. We first study the 
zero-mode background evolution of the spacetime and the inflaton. The background metric 
is 



ds = g^ v dx^dx u = —dt + a (i)cbc 



(2.2) 



where a(t) is the scale factor and x is the comoving spatial coordinates. The equations of 
motion are 



H 2 



1 ,1 : 



H 



3M| V 2 
2Ml ' 



o0o + V) 



+ 3H4> + V = 



(2.3) 

(2.4) 
(2.5) 



The first equation determines the Hubble parameter H, which is the expansion rate of the 
universe. The second equation is the continuity condition. The third equation describes the 
evolution of the inflaton. Only two of them are independent. 

The requirement of having at least O(60) e-fold of inflation imposes some important 
conditions. By definition, to have this amount of inflation, the Hubble parameter cannot 
change much within a Hubble time H" 1 . This gives the first condition 



(2.6) 



We also require that the parameter e does not change much within a Hubble time, 



v = 7h <<c oil) ■ 



(2.7) 



In principle, rj can be close to 0(1) but e kept small. In such a case, e grows exponentially 
with efolds and the inflation period tends to be shorter. More importantly, such a case will 
not generate a scale-invariant spectrum, as we will see shortly, thus cannot be responsible for 
the CMB. The two conditions (I2.6P and (12. 7p are called the slow-roll conditions. Using the 
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background equations of motion, we can see that the slow-roll conditions impose restrictions 
on the rolling velocity of the inflaton. The first condition ( \2.6\\ implies that 



So the energy driving the inflation on the right hand side of (I2.3P is dominated by the 
potential. Adding the second condition (I2.7P further implies that 

°" -e+J<0(l). (2.9) 



So the first term O i n (12. 5 p is negligible and the evolution of the zero-mode inflaton is 
determined by the attractor solution 

3#0 O + V = . (2.10) 

Using (I2.10p . the slow-roll conditions can also be written in a form that restricts the shape 
of the potential, 

M 2 fV'\ 2 V" 

They are related to e and rj by 

e = e y , 7] = -2r] V + Ae v . (2.12) 



So the shape of the potential has to be rather flat relative to its height. We emphasize that, 
although in this example several definitions of the slow-roll conditions are all equivalent, the 
definition (12. 6p and (12. 7p are more general. In other cases that we will encounter later in this 
review, these two conditions are still necessary to ensure a prolonged inflation and generate 
a scale-invariant spectrum, but the others no longer have to be satisfied. For example, the 
shape of potential can be steeper, or the inflationary energy can be dominated by the kinetic 
energy. 

Now let us study the perturbations. To keep things simple but main points illustrated, 
in this section, we will ignore the perturbations in the gravity sector and only perturb the 
inflaton, 



0(x,t) = o (t) + 50(x,t) . (2.13) 

We also ignore terms suppressed by the slow-roll parameters, which we often denote collec- 
tively as 0(e). For example, the mass of the inflaton is V" ~ 0(e)H 2 and will be ignored. 
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The quadratic Lagrangian for the perturbation theory is simply 

L = I d 3 x 



(2.14) 



a -H 2 - a -{dM? 



and the equation of motion follows, 



k 2 

5(p(k, t) + 3H5^(k, t) + -^(k, t) = . (2.15) 



a 2 



where we have written it in the comoving momentum space, 

5(f)(k,t) = J d 3 a^(x,t)e ikx . (2.16) 

The solution to the differential equation (I2.15p . u(k, t), is called the mode function. It is not 
difficult to check that 

a 3 u(k, t)u*(k, t) — c.c. = t-independent const. . (2-17) 

To quantize the perturbations according to the canonical commutation relations between 8(f) 
and its conjugate momentum 5ti = dL/d5<p, 

[50(x, t), 5vr(y, t)] = itf(x - y) , 

[<50(x, t), 50(y, t)] = Q, [5tt(x, t), 5n(y, t)] = 0, (2.18) 



we decompose 



5<p = w(k,t)a k + u*(-k,t)at k , (2.19) 
5tt = a 3 tt(k, t)a k + a 3 w*(-k, t)a^ k , (2.20) 



with the commutation relations 

[a p ,aL q ] = (27r) 3 5 3 (p + q), 

[a p ,a_ q ]=0, [a p ,a f _ q ] = 0. (2.21) 

One can check that the commutation relations (I2.18P and (12.211) are equivalent because of 
(12.17p . given that the constant on the right hand side of (I2.17P is specified to be i. This gives 
the normalization condition for the mode function. 

We now write down the mode function explicitly by solving (I2.15p . 

u(k, r) = C + -^=(l + ikr)e~^ + CL-|=(1 - ikr)e ikT , (2.22) 



7 



where we have used the conformal time r defined as dt = adr. The infinite past corresponds 
to r — > — oo and the infinite future r — )■ 0. We also used the relation r = — l/i/a+(9(e). This 
mode function is a superposition of two linearly independent solutions with the normalization 
condition 

|C + | 2 - |C_| 2 = 1 (2.23) 

followed from (12.171) . Consider the limit in which the mode is well within the horizon, i.e. its 
wavelength a/k much shorter than the Hubble length 1/H, and consider a time period 
much shorter than a Hubble time. In these limits, the mode effectively feels the Minkowski 
spacetime, and the first component in ( 12. 22ft with the positive frequency asymptotes to 
the vacuum mode of the Minkowski spacetime as we can see from (I2.23[) . We choose this 
component as our vacuum choice, and it is usually called the Bunch-Davies state. The 
annihilation operator a p annihilates the corresponding Bunch-Davies vacuum, a p |0) = 0. 
The mode function 

u(k, r) = ^=(1 + ikr)e- ikT (2.24) 

has the following important properties. It is oscillatory within the horizon k\r\ 3> 1- As it 
gets stretched out of the horizon k\r\ <C 1, the amplitude becomes a constant and frozen. 
Physically this means that, if we look at different comoving patches of the universe that 
have the superhorizon size, and ignore the shorter wavelength fluctuations, they all evolve 
classically but with different 5(f). This difference makes them arrive at 0/, the location of the 
end of inflation, at different times. This space-dependent time difference 5t « <50/<A) leads 
to the space-dependent inflationary e-fold difference 

C^HSt^ hJ- . (2.25) 

Again we ignore terms that are suppressed by the slow-roll parameters. This e-fold dif- 
ference is the conserved quantity after the mode exits the horizon, and remains so until 
the mode re-enters the horizon sometime after the Big Bang. It is the physical quantity 
that we can measure, for example, by measuring the temperature anisotropy in the CMB, 
( « —5AT/T [39] . The information about the primordial inflation is then encoded in the 
statistical properties of this variable. So we would like to calculate the correlation func- 
tions of this quantity. Using ( 12.251) . ( 12.191) . ( 12.241) and ( 12. 81) . we get the following two-point 
function, 

(C 2 ) = (0|C(k 1 ,0)C(k 2 ,0)|0) = ||(27r) 5 5(k 1 + k 2 ) , (2.26) 
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where is defined as the power spectrum and in this case it is 



p < = Jki ■ (2 - 27) 

The spectrum index is defined to be 

d\nP c 

id = - 2£ -"' (2 ' 28) 

where the relation dink = Hdt is used. If n s = 1, the spectrum is scale invariant. The 
current data from CMB tells us n s = 0.963 ± 0.012 [1J. So as we have mentioned, this 
requires a small rj, which is also a value that tends to give more e-folds of inflation. 

If this were the end of story, all the primordial density perturbations would be determined 
by this two-point function and they are Gaussian. The rest of the review will be devoted to 
making the above procedure rigorous and to the calculations of higher order non-Gaussian 
correlation functions in this and various other models. 



3 In-in formalism and correlation functions 

In this section, we review the in-in formalism and the related techniques that are used to 
calculate the correlation functions in time-dependent background. The main procedure is 
summarized in the last subsection. 

3.1 In-in formalism 

We start with the in-in formalism |40fH5] . following Weinberg's presentation [35] . 

We are interested in the correlation functions of superhorizon primordial perturbations 
generated during inflation. So our goal is to calculate the expectation value of an operator 
Q, which is a product in terms of field perturbations 5<ft a and 5ir a , at the end of inflation. 
The subscript a labels different fields. In inflation models, these fields are, for example, 
the fluctuations of the scalars and metric and their conjugate momenta. In the Heisenberg 
picture, 

(Q) = (Q\Q(t)\Q) , (3.1) 

where t is the end of inflation, \Q) is the vacuum state for this interacting theory at the far 
past t Q . 

We start by looking at how the time-dependence in Q(t) is generated. 
The Hamiltonian of the system 

H [<f>(t), 7r(t)] = / d 3 xn [0 o (x, t), 7T a (x, t)) (3.2) 
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is a functional of the fields a (x, t) and their conjugate momenta tt (x, t) at a fixed time 
t. On the left hand side of (13. 2h we have suppressed the variable x and index a which are 
integrated or summed over. The a (x, t) and 7r a (x, t) satisfy the canonical commutation 
relations 



[0 a (x, t), vr fe (y, t)] = i5 ab 5 3 (x - y) , 
[0 a (x,t),0 6 (y,t)] = [7r a (x,t),7r b (y,t)] = , 

and their evolution is generated by H following the equations of motion, 



a ( X) t) = i[ff[0(t),7r(t)],0 o (x,t)] 



7r a (x,t) =i[/r[0(t),7r(t)],7r a (x,t)] 



(3.3) 



(3.4) 



We consider a time-dependent background, a (x, t) and 7r a (x, t) which are c- numbers and 
commute with everything, and the perturbations, 50 a (x, t) and 57r a (x, i), 



a (x, t) = O (X, t) + 50 a (x, t) , 7T a (x, t) = 7T a (x, t) + 57T a (x, t) . 

The background evolution is determined by the classical equations of motion, 



0a(x,t) 



, 7T a (x,t) 



(3.5) 



(3.6) 



The commutation relations (13.31) become those for the perturbations, 



a (x, t), Sir b {y, t)} = i5 ab 5 3 (x - y) , 
[5(f) a (x,t),5(f) b (y,t)] = [5ir a (x,t),5Tt b (y,t)) =0 . 



We expand the Hamiltonian as 

H [<j>(t), 7t(t)] = H 7f(t)] + / rf3x g/(x t) 6M ^ t] + S / ^ 
+ #[<ty(t),MM , 



<9H 



a; 



dir a (x,t) 



(3.7) 



57T a (x,t) 



(3.8) 



where we use if to denote terms of quadratic and higher orders in perturbations. 
Using (13. 6p . (13. 7p and (13.81) . the equations of motion ( 13 .4p become 



«(x,*)=i # [50(t),5vr(t);t],(50 a (x,t) , <57r a (x,t)=i # [<ty(*)» *] , <M X , t) 



(3.9) 



So the evolution of the perturbations, 5<f) a and S7i a , is generated by if. It is straightforward 
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to verify that the solutions for (13. 9p are 

5(/> a (pc,t) = U-\t,t Q )6<f> a (x,t Q )U(t,t ) , 5vr Q (x,t) = f/- 1 (t,to)5vr a (x,to)^(t,to) , (3.10) 
where U satisfies 

jU(t,t ) = -iH[6<f>(t ),67r{t o y,t}U{t,t ) (3.11) 

with the condition at an initial time to being 

U(t ,t ) = l . (3.12) 

To have a systematic scheme to do the perturbation theory, we split H into two parts, 

H [8<f>(t), 5ir(t); t] = H [5(f>(t), 5ir(t); t] + [5<f>(t), 5ir(t); t] . (3.13) 

The Ho is the quadratic kinematic part, which in the perturbation theory will describe 
the leading evolution of fields. Fields whose evolution are generated by H are called the 
interaction picture fields. We add a superscript " J" to label such fields. They satisfy 

6j>l(x,t)=i[H [5<p\t),5^\t)-t\ ,^(x,t)] , 6*i(x,t) =i [H [S^^M^t] MM] 

(3.14) 

The solutions are 

^(x,t) = Uo\t,to)6<f> a (x,t )U (t,t ) , &irl(x,t) = U^^^Tta^^Uoit,^) , (3.15) 
where U satisfies 

jU (t, t ) = -iH [6<f>(to), 8n(to)', t) U (t, t ) (3.16) 

with 

Uo(t ,t ) = l. (3.17) 

So the idea is to encode the leading kinematic evolution in terms of the interaction picture 
fields, and calculate the effects of the interaction through the series expansion in terms of 
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powers of Hj. To do this, we rewrite ( 13.1 j) as 



(n|Q[^ a (x,t),Mx,*)]|n) 

(filZ/- 1 ^,^) g[50 a (x,t O ),57T a (x,t O )] c/(t,t )|fi) 

(n\F-\t, t )Uo% t ) Q [«ty a (x, t ), 5vr a (x, t )] E/ (*, *o)^(*, t )\Sl) 

(n\F-\t,to) Q [<^(x,t),5^(x,t)] F(t,t )|Q) , 



where 



F(Mo) = Uo%t )U(t,t Q ) . 
Using flXTTj) . flXTol) and (I3TT31 . we get 

j t F{t,t ) = - iU Q 1 (t,t )H I [S(P(t ),57r(t y i t}Uo(t,t ) F{t,t ) 

= -iH I [5<p I (t),57c I (t) ] t] F(t,t ) 
= -tHj(t)F(t,t ) , 

with 

F(t ,to) = 1 • 

The solution to (13.201) and (I3.2ip can be written in the following way, 



F(t,t ) = Texp ( -i j H T (t)dt 



(3.18) 



(3.19) 



(3.20) 



(3.21; 



(3.22) 



where the operator T means that, in each term in the Taylor series expansion of the expo- 
nential, the time variables have to be time-ordered. The operator T will be used to mean 
the reversed time-ordering. 

In summary, the expectation value (13. ip is 



(Q) =(Q\F-\t,t )Q I (t)F(t,t )\Q) 



=(Q\ 
Notice that in 



Texp jf" H^db^j Q 1 '(t) T exp ^-i jf Hj(t)dtj 



H!{t) = HT^it)^^ 1 ^] , 



\Q) 



(3.23) 



(3.24) 
(3.25) 



all the field perturbations are in the interaction picture. 
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The perturbation theory is also often done in terms of the Lagrangian formalism. In 
the following, we show that they are equivalent. In the above, we perform perturbations on 
the Hamiltonian, and define 5n a by perturbing i\ a = dL/d<ft a , (here we use d to denote the 
functional derivative,) 

dL ■ dL - - 

S^a = —r(<l>a, <Pa) — (0a, 0a) • (3.26) 

#0a 8(j) a 

The Hamiltonian H is defined by (13. 8p . So using the definition 

r dL 

H = / d 3 x^<p a - L , (3.27) 



together with the classical equations of motions (13. 6p and ix a = dL/d(f> a , the definition (13. 8p 
for H becomes 

/dL f dL - - - - 

d 3 x- r {<j) a , a )50 a + / d 3 x- r {0 a , a )<50 a - L(0 a , <j> a ) + L(0 a , <j> a ) . (3.28) 
d(j) a J dcj) a 

If we perturb the Lagrangian directly, we keep the part of the Lagrangian that is quadratic 
and higher in perturbations S(j) a and 5<j) a , 

/dL - - C dL - - 

d 3 X—^((j) a ,(j) a )8(f) a - I d 3 X^((j) a ,(j) a )5<j) a . 
d<Pa J d(j) a 

(3.29) 

The Sir a is defined directly as 

dL dL ■ dL - - 

6rCa = ^771 n = 7r^ ffi ' ~{<i>a, 0a) , (3.30) 

d(5(f) a ) d<p a d(j) a 

where in the second step Eq. (I3.29P has been used. So these two definitions of 5ir a are 
equivalent. The Hamiltonian H is defined through L, 

r qJ 

H = I d 3 x — ^84> a - L . (3.31) 



Again, using (I3.29j) . we can see that the two definitions of H are equivalent. 
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3.2 Mode functions and vacuum 



The Hamiltonian H in the above formalism is typically chosen to be the quadratic kinematic 
terms for field perturbations 5<f) a without mixing, 



H = J d 3 x 



(3.32) 



So they describe free fields propagating in the time-dependent background. The A, B and 
C are some time- dependent background fields, and they are all positive. The solutions to 
the equations of motion (I3.14p in momentum space, u a (k,t), are called the mode functions, 
where k denotes the comoving momentum. They satisfy the Wronskian condition 

Au a (k, t)ii*(k, t) — c.c. = i , (no sum over a) . (3.33) 

Note that we have specified the time-independent constant on the right hand side of f !3.33|) 
to be % for the same reason that we see in Sec. |IZJ Namely, we decompose as 

<J#(k, t) = u a (k, t)a a {k) + <(-k, *)a+(-k) , (3.34) 

where the annihilation and creation operators satisfy the following relations, 

[a (k),aj(-p)] = (27i) 3 5 ab S 3 (k + p) , 

[a a (k), a 6 (-p)] = , [4(k), aj(-p)] = . (3.35) 

These commutation relations are equivalent to (13. 7p because of (I3.33p . but the constant needs 
to be i. This gives the normalization condition for the mode functions. 

Being the solutions of the second order differential equation, generally the mode function 
is a linear superposition of two independent solutions. So we need to specify the initial 
condition. For inflation models, as long as the field theory applies, one can always find an 
early time at which the physical momentum of the mode is much larger than the Hubble 
parameter and study a time interval much less than a Hubble time. Under these conditions, 
the equations of motion approach to those in the Minkowski limit, in which the mode function 
is a linear superposition of two independent plane waves, one with positive frequency and 
another negative. The ground state in the Minkowski spacetime is the positive one. The 
mode function which approaches this positive frequency state in the Minkowski limit is 
called the Bunch-Davies state. In physical coordinates, this limit is proportional to e~ lkpht , 
(for k p h ^> m), where k p h is the physical momentum. In terms of the conformal time 
t = J dt/a{t) and the comoving momentum coordinate k = k p \Ja{t) which we often use, 
this limit is proportional to e~ tkr . We have seen an example in Sec. |2] and will see more 
similar examples later with different A, B and C. The corresponding vacuum |0) is the 
Bunch-Davies vacuum and annihilated by a a (k) defined in f)3.34p . a a (k)|0) = 0. 
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We also would like to write the vacuum of the interacting theory (I3.23P in terms of 
the vacuum of the free theory |0) defined above. Unlike the scattering theory where the 
vacuum of the free theory is generally different from the vacuum of the interaction theory, 
the process that we are studying here do not generate any non-trivial vacuum fluctuations 
through interactions. This is a direct consequence of the identity 

F~ X F = 1 . (3.36) 

So we can replace \Q) in ( I3.23|) with the Bunch-Davies vacuum |0) that we have specified 
above. 

The integrand Hj(t) in (13. 22ft is highly oscillatory in the infinite past due to the behavior 
of the mode function oc e~ lkr . Their contribution to the integral is averaged out. For the 
Bunch-Davies vacuum, this regulation can be achieved by introducing a small tilt to the 
integration contour r — > — oo(l +ie) or performing a Wick rotation r — > it. The imaginary 
component turns the oscillatory behavior into exponentially decay, making the integral well 
defined. 

3.3 Contractions 

When evaluating (I3.23p . one encounters (anti-)time-ordered integrals, of which the integrands 
are products of fields, such as <50„ and 57rf , or <50„ and 50^, sandwiched between the vacua. 
In contrast to the Minkowski space, in the inflationary background, we do not have a simple 
analogous Feynman propagator which takes care of the time-ordering. Therefore we will 
just evaluate the integrands, but leave the complication of the time-ordering to the final 
integration. 

To evaluate the integrand, one can shift around the orders of fields in that product, 
following the rules of the commutation relations. A contraction is defined to be a non-zero 
commutator between the following components of two fields, [5(p^, 50^] , where 50+ and 5cf)~ 
denote the first and second term on the right hand side of (I3.34p . respectively. After normal 
ordering, namely moving annihilation operators to the right-most and creation operators to 
the left-most so that they give zeros hitting the vacuum, it is not difficult to convince oneself 
that the only terms left are those with all fields contracted. Feynman diagrams can be used 
to keep track of what kind of contractions are necessary. 

In the following we demonstrate this using a simple example. We consider a field 8(f) 1 
and quantize it as usual, 

5(f) 1 (k, t) = 5<f) + + 5(f)' = u(k, t)a k + u*(-k, t)a f _ k . (3.37) 

So a contraction between the two terms, 5(f>(k, t') on the left and 5(p(p,t") on the right, is 
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Figure 1: An example of Feynman diagram. 



defined to be 



[8<f>- 



k,0,ty-(p,Ol 



M (k,t'K(-p,0(2vr) 3 5 3 (k + p) . 



(3.38) 



For example, we want to compute a contribution to the four-point function (<50 4 ) from a 
tree-diagram containing two three-point interactions of the following form, 



These two iJ 7 's come from expanding F _1 or F in ( I3.23p . The corresponding Feynman 
diagram is Fig. [TJ 

In Fig. [TJ the two cubic vertices each represent the three-point interaction (I3.39p . Each 
line represents a contraction. The four outgoing legs connect to the four 5(f>(pi,t) (i = 
1,2,3,4) in (50 4 ). The following is a term from the perturbative series expansion of (I3.23p . 
We demonstrate in the following one set of contractions represented by the diagram in Fig. [U 



= [54>+ (pi , *'), &<F (ki , t)] [64>+ ( P2 , t'), 8<j>- (k a , f)] N + (k 3 , t), 8j>- (qi , *")] [H + (k 4 , *) , 6jr (q 2 , t")) 
[5^ + (p 3 ,t'),S4>-(q 3 ,t")] ■ 



Note that all terms are contracted. The result can be further evaluated using (I3.38p . After 
integration over momenta indicated in ( I3.39p . the final momentum conservation will always 
manifest itself as {2ir) 3 ^j(kj). There are other sets of contractions represented by the 
same diagram for the same term. Namely, there are three ways of picking two of the three 
Pi's (qi's), so we have a symmetry factor 9; also, there are 24 permutations of the four 
kj's. We need to sum over all these possibilities. We also need to sum over all possible 
terms containing two if 7 's in the perturbative series, which are not listed here, with their 
corresponding time-ordered integral structure. 




(3.39) 



6^ (pi , t')S^ ( P2 , TjsJ 1 ( P3 , (ki , (k 2 , t)6<? (ka , t)8<f (k 4 , t)5^ ( qi ,t")8^ ((i2,t")d^ (q 3 ,t") 
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3.4 Three forms 



Now we deal with the time-ordered integrals in the series expansion. There are two ways to 
expand (13.231) . 

In the first form, we simply expand the exponential in (I3.22p . For example, for an even 
n, the nth order term is 

ft /"fa fin/2-l ft ftl ftn/2-1 

i n (-l) n/2 dh dt 2 --- dt n/2 dh dt 2 --- dt n/2 

Jt() Jto J U) Jto Jto Jto 

x (iZ/(t„/2) ■ ■ • HtfjQiWHjfa) ■ ■ ■ Hj(t n/2 )} 

n/2 



f - f - ft-n/2-l-m _ ft ftl ftn/2-l + m 

+2Re^>(-l) m+n/2 / dh / dt 2 --- di n /2-m dh dh-- dt n/2+ri 

m= l Jto Jto Jto Jto Jto Jto 

x (Hi(t n/2 „ m ) ■ ■ ■ Hj^Qi^Hjih) ■ ■ ■ Hi(t n/2+m )) . (3.40) 

Each term in the above summation contains two factors of multiple integrals, one from F^ 1 
and another from F. Each multiple integral is time-ordered or anti-time-ordered, but there 
is no time-ordering between the two. We call this representation the factorized form. 

In the second form, we rearrange the factorized form so that all the time variables are 
time-ordered, and all the integrands are under a unique integral. The nth order term in this 
form is [45] 

ft fh ftn-l 

i n dh dt 2 --- dtndHjiQ^Hrit^),--- Jif^),^)] •••]]) . (3.41) 

Jto Jto Jto 

We call this representation the commutator form. 

Each representation has its computational advantages and disadvantages. 

The factorized form is most convenient to achieve the UV (t{ — > to) convergence. As 
mentioned, after we tilt or rotate the integration contour into the positive imaginary plane 
for the left integral, and negative imaginary plane for the right integral, all the oscillatory 
behavior in the UV becomes well-behaved exponential decay. However this form is not 
always convenient to deal with the IR (tj — > t) behavior. For cases where the correlation 
functions have some non-trivial evolution after modes exit the horizon, as typically happens 
for inflation models with multiple fields, the convergence in the IR is slow. Cancellation of 
spurious leading contributions from different terms in the sum (I3.40p can be very implicit 
in this representation, and could easily lead to wrong leading order results in analytical 
estimation or numerical evaluation. 

The commutator form is most convenient to get the correct leading order behavior in 
the IR. The mutual cancellation between different terms are made explicit in terms of the 
nested commutators, before the multiple integral is performed. However, such a regrouping 
of integrands makes the UV convergence very implicit. Recall that the contour deformation 
is made to damp the oscillatory behavior in the infinite past. In the commutator form, for 
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any individual term in the integrand, we can still generically choose a unique convergence 
direction in terms of contour deformation. Although the directions are different for different 
terms, they can be separately chosen for each of them. But now the problem is, if these 
different terms have to be grouped as in the nested commutator so that the IR cancellation 
is explicit, the two directions get mixed. Hence the explicit IR cancellation is incompatible 
with the explicit UV convergence in this case. 

To take advantage of both forms, we introduce a cutoff t c and write the IR part of the 
in-in formalism in terms of the commutator form and the UV part in terms of the factorized 
form [4"6] . 

n ft fU-l ftc ft„-l 

/ dt\ • ■ ■ I dti {commutator form} / dti + \ ■ • ■ I dt n {factorized form} . 

J tc J tc J —CO J —CO 

(3.42) 

This representation is called the mixed form. This form is particularly efficient in numerical 
computations when combined with the Wick-rotations in the UV. 

We will not always encounter all these subtleties in every model, but there does exist 
such interesting examples, as we will see in Sec. 17.11 



3.5 Summary 

To end this section, we summarize the procedure that we need to go through to calculate 
the correlation functions in the in-in formalism. 

Starting with the Lagrangian L[<f)(t), <j>(t)), we perturb it around the homogenous solutions 
4> a and (j) a of the classical equations of motion, 



(/> a (x, t) = 4> a {t) + <50 a (x, t) , a (x, t) = 4> a (t) + <50 a (x, t) . 



(3.43) 



Keep the part of the Lagrangian that is quadratic and higher in perturbations and denote 
it as L. Define the conjugate momentum densities as 57t a = dL/d(5<j) a ). We can also 
equivalently expand the Hamiltonian H[<f>(t),Tr(t)] by perturbing <f> a (pc,t) and 7r a (x, t). 

Work out the Hamiltonian in terms of 5<p a and 5n a , and separate them into the quadratic 
kinematic part H , which describes the free fields in the time-dependent background, and the 
interaction part Hi. Relabel <50 a 's and £7r a 's in the Hamiltonian density as the interaction 
picture fields, Scp^s and Sn^s. These latter variables satisfy the equations of motion followed 
from the H Q . We quantize £0^ and <57r^ in terms of the annihilation and creation operators 
as in (13.341) and (13.351) . The mode functions u a (k, t) are solutions of the equations of motion 
from Hq, normalized according to the Wronskian conditions f!3.33|) and specified by an initial 
condition such as the Bunch-Davies vacuum. The correlation function for Q(t) is given by 



(Q(t)) = (o| 



Texp 



t J 



Q\t) 



Texp ( -% I H^dt 

to 



10} 



(3.44) 
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where Q(t) is a product in terms of <50 7 (x, t) and 57r 7 (x, t). If we want to work with 5^ and 
6^ instead of 5^ and 5tc^, we replace dn^ with using the relation 5^ = dH /d(S7i^). 

Choose appropriate forms in Sec. 13.41 and series-expand the integrand in powers of Hj 
to the desired orders. Perform contractions defined in Sec. 13.31 for each term in this expan- 
sion. Each term gives a non-zero contribution only when all fields are contracted. Draw 
Feynman diagrams that represent the correlation functions, and use them as a guidance to 
do contractions. Finally sum over all possible contractions and perform the time-ordered 
integrations. 



4 A no-go theorem 

Simplest inflation models generate negligible amount of non-Gaussianities that are well below 
our current experimental abilities [UJHH]. By simplest, we mean 

• single scalar field inflation 

• with canonical kinetic term 

• always slow-rolls 

• in Bunch-Davies vacuum 

• in Einstein gravity. 

This list is extracted based on Maldacena's computation of three-point functions in an ex- 
plicit slow- roll model jlT]. We now review this proof. The notations here follow Ref. [4~9|l50] 
and will be consistently used later in this review. 

The Lagrangian for the single scalar field inflation with canonical kinetic term is the 
following, 



S = d A xy/^g 



(4.1) 



where is the inflaton field, X = —^g^ u d^(pd u (f) is the canonical kinetic term and V is the 
slow-roll potential. The first term is the Einstein gravity and Mp = (87rG) _1 / 2 is the reduced 
Planck mass. For convenience we will set the reduced Planck mass Mp = 1. The signature 
of the metric is (—1, 1, 1, 1). 

The inflaton starts near the top of the potential and slowly rolls down. As we have 
reviewed in Sec. [2j to ensure that the inflation lasts for at least O(60) efolds, the potential 
is required to be flat so that the slow-roll parameters (12.111) are both much less than one 
most of the time. The energy of the universe is dominated by the potential energy, and the 
inflaton follows the slow-roll attractor solution (I2.10p . Also as discussed in Sec. [2j we will 
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use the following more general slow-roll parameters, 



e -~W' * = 7h- (42) 

To study the perturbation theory, it is convenient to use the ADM formalism, in which 
the metric takes the form 

ds 2 = -N 2 dt 2 + h ij (dx i + N i dt)(dx j + N j dt) . (4.3) 

The action becomes 

S = ^J dtdx 3 VhN(R^ + 2X - 2V) + ^ j dtdx^VkN-^E^ - E 2 ) , (4.4) 

where the index of N % can be lowered by the 3d metric hij and is the 3d Ricci scalar 
constructed from h^. The definition of E^ and E are 

Eij = lihij-ViNj-VjNi) , 
E = Eijh ij . (4.5) 

In the ADM formalism, the variables N and N l are Lagrangian multipliers whose equations 
of motion are easy to solve. In single field inflation, we have only one physical scalar pertur- 
bation [51] . We choose the uniform inflaton gauge (also called the comoving gauge) in which 
the scalar perturbation ( appears in the three dimensional metric hij in the following form, 

h^ = a 2 e 2C, 5ij , (4.6) 

and the inflaton fluctuation 5(f) vanishes. The a(t) is the homogeneous scale factor of the 
universe, so ( is a space-dependent rescaling factor. In this review we do not consider the 
tensor perturbations. 

We plug (14 .3p and (14.61) into the action (14.41) and solve the constraint equations for the 
Lagrangian multipliers N and N l . We then plug them back to the action and expand up 
to the cubic order in ( in order to calculate the three-point functions. To do this, in the 
ADM formalism, it is enough to solve N and N l to the first order in (. This is because their 
third order perturbations will multiply the zeroth order constraint equation which vanishes, 
and their second order perturbations will multiply the first order constraint equation which 
again vanishes. After some lengthy algebra, we obtain the following expansions, 



So = / dtd 3 x 



a 3 e( 2 - ae(d() 2 , (4.7) 
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where 



dtd 3 x 



a s e\( 2 + ae 2 ((d(y 2 -2ae((d0(dx) 



+ /(C) 



5L 



(4- 



X 

6L 

/(C) 



a 2 ed- 2 ( , 

dd 2 x 
dt 



2a 



+ Hd\ - ed 2 ( 
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—C 2 + terms with derivatives on C 
4 s s 



(4.9) 
(4.10) 

(4.11) 



Here d~ 2 is the inverse Laplacian and SL/S(\i is the variation of the quadratic action with 
respect to the perturbation (. We now can follow Sec. [3] and proceed to calculate the 
correlation functions. For simplicity, we will always neglect the superscript "J" on various 
interaction picture fields. 

We restrict to the case where the slow-roll parameters are always small and featureless. 
We first look at the quadratic action. In this case, we can analytically solve the equation of 
motion followed from (14. 7p in terms of the Fourier mode of (, 



Uk 



d 3 x((t, x)e" 



ik-x 



(4.12) 



and get the mode function 



u k = u(k,r) 



iH 



(1 + ikr)e 



-iki 



(4.13) 



where r = J dt/a ~ —(aH)~ l is the conformal time. The normalization is determined by 
the Wronskian condition (I3.33p . We have chosen the Bunch-Davies vacuum by imposing the 
condition that the mode function approaches the vacuum state of the Minkowski spacetime 
in the short wavelength limit k/a 3> ^/H, 



u k ->■ 



Ht 



-ikr 



'Aek 



(4.14) 



The dynamical behavior of ( that has been explained around Eq. (I2.24p and ( I2.25|) is made 
precise here. In particular, ( is exactly massless without dropping any C(e) suppressed 
terms. In addition, from (14. 6p . we can see that, for superhorizon modes, the only effect of ( 
is to provide a homogeneous spatial rescaling. And ( is the only scalar perturbation. So the 
fact that ( is frozen after horizon exit will not be changed by higher order terms. 
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If we choose the spatially flat gauge, we make ( disappear and the scalar in this pertur- 
bation theory becomes the perturbation of <fi. The relation between ( and 5(f) in (I2.25f) (with 
0(e) corrections) is thus a gauge transformation through a space-dependent time-shift. 

We quantize the field as 



C(k, r) = w fc a k + u* k a\ , (4.15) 

with the canonical commutation relation [a^, a^,] = (27r) 3 <5 3 (k — k'). We can easily compute 
the two-point function at the tree level, 



(C(k 1 )C(k 2 )) = ^(2 7 r) 5 5 3 (k 1 + k 2 ) , (4.16) 

where 



p < = i§c- IU71 



5 n^e 



Since ( remains constant after it exits the horizon, the H and e are both evaluated near the 
horizon exit. 

We next look at the cubic action. For single field models, Hj^ = — L 3 . Keeping in mind 
that x is proportional to e, one can see that the first line of f 14 . 8 [) is proportional to e 2 . For 
the featureless potential, fj = 0(e 2 ), where e collectively denotes either e or rj. So the second 
line of (14. 8p is proportional to e 3 , and negligible. The third line can be absorbed by a field 
redefinition ( — > ( n + f(( n ). The only term in f(( n ) that will contribute to the correlation 
function is written out explicitly in (14. lip . All the others involve derivatives of ( so vanish 
outside the horizon. Thus this redefinition eventually introduces an extra term 

(C(k 1 )C(k 2 )C(k 3 )) = (Cn(k!)Cn(k 2 )Cn(k 3 )) 

+ |((C^(k 1 )Cn(k 2 )C„(k 3 )) + 2 perm.) + 0fo 2 (p£ ) 3 ) . (4.18) 

According to (I3.44p . we expand the exponential to the first order in H I 3 to get the leading 
result, 

(O = -<<0| / < ^[Cn(k 1 )Cn(k 2 )C„(k 3 ),^ 3 ]|0) . (4.19) 

To estimate the order of magnitude of the bispectrum, we only need to keep track of 
the factors of H and e. For example, from the first term in (14. 8p . we have J dtH^{t) D 
— J dx 3 d,Ta 2 e 2 ((' 2 , where we used the conformal time r and the prime denotes the derivative 
to r. Using a oc H~ x , ( oc (' oc H/yfe, we see that this three-point vertex contributes 
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cx Hy/e. Together with the three external legs ( 3 and the definition oc H 2 /e, we get 



<0 = -'' / ^/(c 5 . ///..(/)) v— vO(o/; 



(4.20) 



Similar results can be obtained for the other two terms in the first line of ( 14. 8ft . As we will 
define more carefully later, the size of the three-point function is conventionally characterized 
by the number f NL , which is defined as (( 3 ) ~ fmPf- So the contribution from the first line 
of (14.81) is fNL = 0(e). The extra term due to the redefinition (14.181) contributes f^L = 0(rj). 
This completes the order-of-magnitude estimate. To get the full non-Gaussian profile, we 
need to compute the integrals explicitly and get 



(C(k 1 )C(k 2 )C(k 3 )) = (27r) 7 5 3 (k 1 + k 2 + k 3 )P ( 



S 



c n,*r ' 



(4.21) 



where 



S 



fc? \ fh \ 8 fhk 2 

— — h 2 perm. ] + I - — ho perm. ] + — I — h 2 perm. 

k 2 k 3 J \k 2 J K \ k 3 



+ ? ( -r^~ + 2 perm. 



8 V k 2 k 



2^3 



(4.22) 



where K = k\ + k 2 + k% and the permutations stand for those among k±, k 2 and k%. 

The slow-roll parameters are of order (9(0.01), so f^L ~ O(0. 01) for these models. Even 
if we start with Gaussian primordial perturbations, non-linear effects in CMB evolution 
will generate /nl ~ 0(1) [37], and a similar number for large scale structures due to the 
non-linear gravitational evolution or the galaxy bias [35]. It seems unlikely that we can 
disentangle all these contaminations and detect such small primordial non-Gaussianities in 
the near future. 



5 Beyond the no-go 
5.1 Inflation model building 

The following are two examples of slow-roll potentials in the simplest inflation models that 
we studied in Sec. HJ 

V - im 2 2 , (5.1) 
^m 2 2 . (5.2) 

The first type ( 15. ip belongs to the small field inflation models. The slow-roll conditions 
(12. lip require the potential to be flat enough relative to its height, i.e. the mass of the 



^small 
^large 
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inflaton should satisfy m <C H . The second type (j5.2p belongs to the large field inflation 
models. The potential also needs to be flat relative to its height, but here one achieves this 
by making the field range very large, typically ^ Mp. The other conditions that we 
listed in Sec. H] should also be satisfied by these models. These are the classic examples, 
which exhibit algebraic simplicities and illustrate many essential features of inflation. 

However, when it comes to the more realistic model building in a UV complete setup, such 
as in supergravity and string theory, situations get much more complicated. For example, it 
is natural that we encounter multiple light and heavy fields, and the potentials for them form 
a complex landscape. These multiple fields live in an internal space, whose structure can be 
very sophisticated. In string theory, some of them manifest themselves as extra dimensions 
and can have intricate geometry and warping. All these elements have to coexist with the 
inflationary background that introduces profound back-reactions. 

Even with varieties of model building ingredients, it has been proven to be very subtle 
to construct an explicit and self-consistent inflation model. Indeed various problems have 
been noticed over the years in the course of the inflation model building. For example, 

• The ^-problem for slow-roll inflation [52J. As we have seen, in order to have 
slow-roll inflation [8j[9], the mass of the inflaton field has to be light enough, m <^ H, 
to maintain a flat potential. However, in the inflationary background, the natural 
mass of a light particle is of order H. This can be seen in many ways, and in some 
ideal situations they are equivalent to each other. For example, one way to see this 
is to consider the coupling between the Ricci scalar and the inflaton, ~ R<p 2 . In the 
inflationary background R ~ H 2 . Unless we have good reasons to set the coefficient of 
this term to be much less than one, it will give inflaton a mass of order H, spoiling the 
inflation. Another way to see this is to note that the effective potential in supergravity 
takes the form V = Voexp(K/Mp) x other terms. Here schematically K ~ 2 + • • • is 
the Kahler potential and its dependence on is normalized as such to give the canonical 
kinetic term for 0. So the first term in the expansion of V is of order Vo0 2 /M| ~ H 2 <f) 2 
and model independent. Therefore, either symmetry needs to be imposed or other 
tuning contributions introduced to solve this 77-problem. 

• The /i-problem for DBI inflation [S3]- DBI inflation [53] is invented to generate 
inflation by a different mechanism. It makes use of the warped space in the internal 
field space [55JI56]. These warped space impose speed-limits for the scalar field, so 
even if the potential is steep, the inflaton is not allowed to roll down the potential very 
quickly. A canonical example of warped space is 

ds 2 = h(r) 2 (-dt 2 + a(t) 2 dx 2 ) + h(r)- 2 dr 2 , (5.3) 

where r is the extra-dimension (or internal space), h(r) = r/R is the warp factor and R 
is the length scale of the warped space. The position of a 3 + 1 dimensional brane in the 
r-coordinate is the inflaton. So the inflaton velocity is limited by the speed-limit in the 
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r-direction, h 2 . In order to provide a speed-limit that is small enough for inflation, the 
warp factor has to be small enough, h <C HR. However one of the Einstein equations 
with the metric (15.31) takes the following form, 

(dh/dr) 2 — H 2 h~ 2 = — + other source terms , (5.4) 

R 2 

where the second term on the left hand side is due to the back-reaction of the infla- 
tionary spacetime. It is easy to see that the naive h = r/R should be modified for 
h < HR, precisely where the inflation is supposed to happen. Without contributions 
from other source terms, such a deformed geometry closes up too quickly and leads 
to an unacceptable probe-brane back-reaction if we demand the inflaton still follow 
the speed-limit. Therefore, either symmetry, or tuning using other source terms from 
the right hand side of (15 .4p . is necessary to solve this /i-problem. The ^-problem and 
/i-problem are closely related in an AdS/CFT setup. 



The field range bound [57| [58] . Large field inflation models require the field range 
to be much larger than Mp. In supergravity and string theory, starting from a ten- 
dimensional theory with 10-dim Planck mass M 10 , the 4-dim Planck mass Mp is ob- 
tained by integrating out the six extra-dimensions, 

Mf 10) J d 6 yd 4 x^G^R m 
D M ( 8 10) V (6) J d 4 x^^R w = M 2 J d 4 Xy/^R {4) , (5.5) 

where we use L and V( 6 ) ~ L 6 to denote the size and volume of the extra-dimensions, 
respectively. The field range A<p often appears as the distance in the extra dimensions, 
Acf) ~ AL ■ M 2 1Q y with the factor M 2 10 ^ being the proportional coefficient. Clearly, 
AL < L. If the field range manifests itself within a warped throat with a length scale 
R, we still require R < L, and so AL < L. Together with M P = M^L 3 , we get 

A0 < M P /(M (10) L) 2 . (5.6) 

We further note that the microscopic length scale L has to be much larger than the 
10-dim Planck length for the field theory to make sense. So M^L ^> 1, and 

the field range A0 in these models is generically sub-Planckian. For example, for a 
warped throat with charge N, (M( 10 )L) 2 > (M( 10 )-R) 2 ~ iV 1 / 2 , we have 

A(j) < Mp/VN . (5.7) 

We have ignored a detailed numerical coefficient appearing on the right hand side of 
(15.7j) . which is model dependent. For example, considering the volume V(e) to be the 
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sum of the throat and a generic bulk volume, it is (9(0.01) [57J; considering an extreme 
case where the throat does not attach to a bulk, it is 0{1) [58J. Notice that, due to the 
dependence of Mp on the volume V^, increasing the volume only makes the bound 
tighter. 

• The variation of potential (59]. Even in cases where there is no fundamental re- 
striction on the excursion of fields, one encounters problems constructing the large field 
inflationary potential. Large field potentials that arise from a fundamental theory take 
the following general from, 



where mf un d represents typical scales in the theory. For field theory descriptions to hold, 
such scales are much less than Mp. For example, mf UIK j can be the higher dimensional 
Planck mass, string mass, or their warped scales. The A n 's are dimensionless couplings 
of order 0(1). Unless some symmetries are present to forbid an infinite number of 
terms in (I5.8p . or a high degree of fine-tuning is assumed, the shape of potential (15. 8p 
varies over a scale of order mf UI1( j <C Mp. This variation is too dramatic for the potential 
to be a successful large field slow-roll potential. 

None of the arguments in the above list is meant to show that the specific type of inflation 
is impossible. In fact, these have been the driving forces for the ingenuity and creativity in 
the field of inflation model building. This list is used to demonstrate some typical examples 
of complexities in reality. Often times, solving one problem will be companied by other 
structures that make the model step beyond the simplest one. So we may want to keep an 
open mind that the algebraic simplicity may not mean the simplicity in Nature. 



Following is a partial list of possibilities that allow us to go beyond the no-go theorem in 
Sec.H 



• Instead of single field inflation, we can consider quasi-single field or multifield inflation 
models (Sec. & E) . 

• Instead of canonical kinetic terms, there are models where the higher derivative kinetic 
terms dominate the dynamics (Sec. 16. ip . 

• Instead of following the attractor solution such as the slow-roll precisely, features can 
be present in the potentials or internal space, that temporarily break the attractor solu- 
tion, or cause small but persistent perturbations on the background evolution (Sec. 16.21 



• Instead of staying in the Bunch-Davies vacuum, other excitations can exist due to, for 
example, boundary conditions or low scales of new physics (Sec. 16. 4p . 



oo 





n=0 
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• Although strong constraints, from experimental results and theoretical consistencies, 
exist on non-Einstein gravities, early universe may provide an opportunity for their 
appearance. We use this category to include a variety of possibilities, such as modified 
gravities, non-commutativity, non-locality and models beyond field theories. 

There are also strong motivations from data analyses for us to search and study different 
large non-Gaussianities. The signal-to-noise ratio in the CMB data is not large enough for 
us to detect primordial non-Gaussianities model-independently. A well established method 
is to start with a theoretical non-Gaussian ansatz, and construct optimal estimators that 
compare theory and data by taking into accounts all momenta configurations. This then 
gives constraints on the parameters characterizing the theoretical ansatz. Therefore, the 
following two important possibilities exist. First, the primordial non-Gaussianities exist in 
data could be missed if we did not start with a right theoretical ansatz. Second, even if a 
non-Gaussian signal were detected with one ansatz, it does not mean that we have found 
the right one. So different well-motivated non-Gaussian templates are needed for clues on 
how corresponding data analyses should be formed. From a different perspective, even if 
the primordial density perturbations were Gaussian, we would still do the similar amount of 
work and reach the conclusion after various well-motivated non-Gaussian forms are properly 
constrained. 

5.2 Shape and running of bispectra 

In this review, we will be mainly interested in the three-point correlation functions of the 
scalar primordial perturbation (. They are also called the bispectra. In this subsection, we 
introduce some simple terminologies that we often encounter in studies of bispectra. 

The three-point function is a function of three momenta, k 1; k 2 and k 3 , which form a 
triangle due to the translational invariance. Assuming also the rotational invariance, we are 
left with three variables, which are their amplitudes, k\, k 2 and k%, satisfying the triangle 
inequalities. The information is encoded in a function S(ki, k 2 , k^) that we define as below, 

i=l 

where is the fiducial power spectrum, and we fix it to be a constant P^ = P((A;„map) = 
6.1 x 1CT 9 , where /c wmap = 0.027Mpc _1 . We have chosen the above definition so that it can 
be uniformly applied to different types of bispectra that we will encounter in this review. In 
literature, different notations have been used. The differences are simple and harmless. For 
example, different functions such as A = k\k 2 k^S or F = S '/ \kik 2 k^) 2 are sometimes defined. 
We choose S since it is dimensionless and, for scale-invariant bispectra, it is invariant under 
a rescaling of all momenta. This quantity is the combination that is used to plot the profiles 
of bispectra in literature any way, despite of different conventions. Also, the precise power 



27 



(a) (b) (c) 

Figure 2: Momentum configurations: (a) equilateral, (b) squeezed, (c) folded. 

spectrum P^ instead of P^ is often used in the definition (15. 9p . Here we absorb the momentum 
dependence of P^ in S. This is because the three-point function is an independent statistic 
relative to the two-point. In cases where both the power spectrum and bispectrum have 
strong scale dependence, it is not convenient if they are defined in an entangled way. 

Under different circumstances, different properties of S are emphasized. The conventions 
involved may not always be precisely consistent with each other, since they are chosen to 
best describe the case at hand. Following are some typical examples. 

The dependence of 5* on k\, k 2 and k 3 is usually split into two kinds. 

One is called the shape of the bispectrum. This refers to the dependence of S on the 
momenta ratio k 2 /ki and k 3 /k\, while fixing the overall momentum scale K = k\ + k 2 + k 3 . 
Several special momentum configurations are shown in Fig. |2j 

Another is called the running of the bispectrum. This refers to the dependence of S on 
the overall momentum scale K = k\ + k 2 + k 3 , while fixing the ratio k 2 jk\ and k 3 /k\. 

For bispectra that are approximately scale- variant, the shape is a more important prop- 
erty J60j[50]. We will encounter such cases in Sec. I6.1[ l7.1l fe l87Ll The amplitude, also called 
the size, of the bispectra is often denoted as /nl by matching 

S(k u k 2 ,k 3 ) kl=k2=ki > %-f NL . (5.10) 

limit 1U 

In this case, Jnl is approximately a constant but can also have a mild running, i.e. a weak 
dependence on the overall momentum K [6UI62]. An index n^c — 1 = d In fNi/d In k is 
introduced to describe this scale dependence. The power spectrum also has a mild running, 
P^ = (k/k ) Tls ~ 1 P ( ^. In this review, when we give explicit forms of S in the approximately 
scale-invariant cases, for simplicity we mostly ignore these mild scale dependence and con- 
centrate on shapes. Shapes of bispectra have been given names according to the overall de- 
pendence of S on momenta. For example, for the equilateral bispectrum, S peaks at the equi- 
lateral triangle limit and vanishes as ~ k 3 jk\ in the squeezed triangle limit (k 3 -C k\ = k 2 ). 
The local bispectrum peaks at the squeezed triangle limit in the form ~ (k 3 /k\)~ l , such 
as the two shape components in (14. 22ft . To visualize the shapes, we often draw 3D plots 
S(l, x 2 , x 3 ), where x 2 and x 3 vary from to 1 and satisfy the triangle inequality x 2 + x 3 > 1. 

There are also cases where the running becomes the most important property, while the 
shape is relatively less important [64^165] . In such cases, the bispectra are mostly functions of 
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K. So f^L defined in (15. 10p has strong scale dependence. Instead, one can choose a constant 
/nl to describe the overall running amplitude. We will encounter such cases in Sec. 16.21 fc 
16.31 In these cases, the shape plot S(l, x 2 , x 3 ) may look nontrivial but this is because it does 
not fix K. 

The above dissection will become less clean for cases where both properties become 
important. 

One thing is clear. The Jnl, that is always used to quantify the level of non-Gaussianities, 
is only sensible with an extra label that specifies, at least qualitatively, the profile of the 
momentum dependence, such as shapes and runnings. 

It is useful to quantify the correlations between different non-Gaussian profiles, because 
as we mentioned in data analyses an ansatz can pick up signals that are not completely 
orthogonal to it. In real data analyses this is performed in the CMB /-space. To have a 
simple but qualitative analogue in the fc-space, we define the inner product of the two shapes 
as 

S ■ S' = / S(ki, k 2 , k 3 )S'(ki, k 2 , k 3 )w(ki, k 2 , k 3 )dk 1 dk 2 dk 3 , (5.11) 
Jv k 

and normalize it to get the shape correlator [601 163] 



S ■ S' 

C{S, S') = ^.^1/2^.^1/2 • ( 5 - 12 ) 
Following Ref. [S3], we choose the weight function to be 

w{h, k 2 , k 3 ) = - — 1 , (5.13) 
ki + k 2 + k 3 

so that the fc-scaling is close to the /-scaling in the data analyses estimator. Later in this 
review, when we use this correlator to estimate the correlations between shapes, we take 
the ratio between the smallest and largest k to be 2/800, close to that in WMAP. A more 
precise correlator should be computed in the /-space in the same way that the estimator is 
constructed. We refer to Ref. [35] for more details. 

In typical data analyses [66H70] , the estimator involves a triple integral of the bispectrum 
over the three momenta fcj. To have practical computational costs, it is necessary that this 
integral can be factorized into a multiplication of three integrals, each involves only an 
individual fcj. This requires the bispectrum to be of the form fi(ki)f 2 (k 2 )f 3 (k 3 ), or a sum 
of such forms. Such a form is called the factorizable form or separable form. The factor 
K~ n may be tolerable since it can be written as (l/r(n)) J °° t n ~ 1 e~ Kt . If the analytical 
result is too complicated, to make contact with experiments we will try to construct simple 
factorizable ansatz or template to capture the main features of the original one. New methods 
that are applicable to non-factorizable bispectrum forms and are more model-independent 
are under active development [71] . 
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6 Single field inflation 



In this section, we relax several restrictions of the no-go theorem on single field inflation mod- 
els and study how large non-Gaussianities can arise. We present the formalisms and compute 
the three-point functions. We emphasize how different physical processes during inflation 
are imprinted as distinctive signatures in non-Gaussianities. Obviously, any mechanism that 
works for single field inflation can be generalized to multi-field inflation models. 



6.1 Equilateral shape: higher derivative kinetic terms 

In this subsection, we study large non-Gaussianities generated by non-canonical kinetic terms 
in general single field inflation models, following Ref. [50J. 

Consider the following action for the general single field inflation 



S = d 4 xy/^g 



(6.1; 



Comparing to (14. ip . we have replace the canonical form X — V with an arbitrary function 
of X = —^g^d^dptp and cf). This is the most general Lorentz- invariant Lagrangian as a 
function of and its first derivative. It is useful to define several quantities that characterize 
the differential properties of P with respect to X [T2"jl4"9]. 

* = P X + 2XP,xx ■ (6 ' 2) 



E = XP tX + 2X 2 P XX = ^ , (6.3) 



A = X 2 P, XX + ^X 3 P, XXX , (6.4) 

where c s is called the sound speed and the subscript "X" denotes the derivative with respect 
to X. The third derivative is enough since we will only study the three-point function here. 

It is a non-trivial question which forms of P will give rise to inflation. The model- 
independent approach we take here is to list the conditions that an inflation model has to 
satisfy, no matter which mechanism is responsible for it. Namely, we generalize the slow-roll 
parameters in (14. 2 p to the following slow- variation parameters 

H e c s 

e = -Jp' v = 7h> s = ^' (6 - 5) 

and require them to be small most of the time during the inflation. The smallness of these 
parameters guarantees the Hubble constant H, the parameter e and the sound speed c s to 
vary slowly in terms of the Hubble time. Similar to the arguments given in the case of 
slow-roll inflation in Sec. [2} these are necessary to ensure a prolonged inflation as well as an 
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approximately scale-invariant power spectrum that we observed in the CMB. 

Following the same procedure that is outlined in Sec. HI we get the quadratic and cubic 
action for the scalar perturbation ( [mHHlEO] ■ The quadratic part is 



So 



did x 



a 3 -( 2 -ae(d(f 



(6.6) 



If the slow- variation parameters are always small and featureless, we can analytically solve 
the equation of motion followed from ( 16. 6ft and get the following mode function, 



iH 



Uk(T) 



[1 + ikc s r)e 



-ikc a 



y/AeCsk 3 

Notice the appearance of c s comparing to (I4.13p . The two-point function is 



(6.7) 



(C(ki)C(k 2 )) = ^(2vr) 5 5 3 (k 1 + k 2 ) 



(6.8) 



with the power spectrum 



H 2 



9>n 2 ec s 



(6.9) 



where the variables are evaluated at the horizon crossing of the corresponding fc-mode. 

To calculate the bispectrum, we look at the cubic action. In the following, we list three 
terms that are most interesting for this subsection, 



5, 



dtd 3 x 



a 3 e 



1 - 



2A 



C 3 + 



3a 3 e 



1 - 



CC 2 - ae 1 



1 



(6.10) 



The full terms can be found in Eq. (4.26)-(4.28) in Ref. 

The order of magnitude contribution from these three terms can be estimated similarly 
as we did in (I4.20p . but now we not only keep factors of H and e, but also factors of c s . Take 
the first term as an example, we write it in terms of the conformal time, 



drH^ir) D / drd i x 



2A 



c 



/3 



(6.11) 



Comparing ( 16 .7p with (I4.13p . we see that there is an extra factor of c s companying r. So we 
estimate dr oc c' 1 and a f» —{Ht)~ 1 oc CsH" 1 . Also, ( oc H/^/ec~ s , but (' oc c s (. Overall, 
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the vertex (16.1 ip contributes 



oc 



H 



'ec s 





2A" 







Multiplying the three external legs ( 3 , and using the definition 



(C 3 ) ~ InlP? 



and P( oc (H/^/ec^) 2 , we get 



The other two terms are similar. A detailed calculation reveals 



(6.12) 



(6.13) 



(6.14) 



(C(k 1 )C(k 2 )C(k 3 )) = (27r) 7 5 3 (k 1 + k 2 + k 3 )(P c ) 



II,* 

^ S Q ~\- S t -\- ~\- S, 



2 

i '"i 



(6.15) 



where we have decomposed the shape of the three-point function into six parts. The first 
two come from the leading order terms that we listed in (I6.10p . 



S r 



- 1 



2A\ Zk^h 
IT J 2K* 



(5-1) (^E^^E^E*. 3 )^ 



(6.16) 
(6.17) 



In terms of /tvl their sizes are 



f x -- 

81 Vc? 



NL 



35 fl 



108 V c 2 s 



- 1 



(6.18) 
(6.19) 



The next four terms come from the subleading terms that we did not list explicitly in (I6.10P 
as well as the subleading contributions from the first two terms. Their orders of magnitude 
are 



NL 



0(4, ^ 



(6.20) 



The detailed profiles can be found in Ref. [50J. 
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Figure 3: Shape of S\ in (I6.16p . 




Figure 4: Shape of S c in ( 16. 17[) . 
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The full results we obtained can be used in different regimes. 



• If we look at the limit, c s < 1 or A/S > 1, the leading order results give two shape 
components, S\ and S c . This result can also be obtained using a simple method of 
considering only the fluctuations in scalar field while neglecting those in gravity [74"1T75] . 
Intuitively, this is because the higher derivative terms are responsible for the generation 
of large non-Gaussianities, and the gravity contribution is expected to be small as we 
saw in Sec. HI Therefore one expands P(X, 0) using 0(x, t) = 0o(£) + <50(x, t). The 
derivatives of P with respect to are ignored because the inflation and scale invariance 
imposes an approximate shift symmetry on P in terms of inflaton 0. We then get two 
terms in the cubic Lagrangian density 



This gives two leading bispectra the same as (I6.16P and (16.171) . The approach that we 
present here gives a rigorous justification to such an method. The subleading order 
component S a may be observable as well. At this limit where the higher derivative 
terms of the inflaton field are dominant, the Lagrangian of the above effective field 
theory are generalized [76J to include, for example, the ghost inflation [77] whose 
Lagrangian cannot be written in a form of P(X, 0). Another two slightly different 
equilateral shapes arise. However it is worth to mention that, generally in single field 
models and Einstein gravity, going beyond P(X, 0) requires adding either terms that 
explicitly break the Lorentz symmetry, or terms with higher time derivatives on 
which cannot be eliminated by partial integration, such as (D0) 2 . Different treatment 
of such terms and discussions on their effects can be found in [78H8U] . 

• If we take the opposite, slow-roll limit, c s — > 1 and A/S — > 0, we recover the two shape 
components S t and S v that we got in Sec. HJ with unobservable size /atl ~ C( e )- 

• We can also look at the intermediate parameter space. In slow-roll inflation models, 
one can also add higher derivative terms J72JHH]. But in order not to spoil the slow-roll 
mechanism, the effect of these terms can only be subdominant. This corresponds to 
c s ~ 1 and A/S < 0(1). Using the full results, we can see that the size of the non- 
Gaussianity is J'nl < Therefore it is important to emphasize that, for the class 
of models we consider here, non-slow-roll inflationary mechanisms, such as the example 
that will be given below, are necessary to generate observable large non-Gaussianities. 

• The other terms that we did not list in (I6.10P (see Ref. [50]) and their canonical limit 
(14 .8p are also useful. These terms are exact for arbitrary values of e, rj and s, so the 
usage of the action is beyond the category of models that we focus on in this subsection. 
As we will see in Sec. I6.2l and Sec. 16.3} it can be applied to the cases of sharp or periodic 
features where these parameters do not always remain small. 



£ 3 = a 




t;Pxx<Po + \p,xxx4>1] S'4> 3 - ^Pxx0o50(V50) 2 . 



(6.21) 
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Figure 5: An ansatz (16. 22ft for the equilateral shape. 
In the rest of this subsection, we focus on the first case. 

In Fig. [3] and HI we draw the shapes of S\ and S c . The two shapes are similar. They both 
peak at the equilateral limit, and behave as S ~ k^/ki in the squeezed limit k% <C k\ = k 2 - 
We call these shapes the equilateral shapes. There are some small differences between S\ 
and S c , e.g. around the folded triangle limit k 2 + k% = k\. A factorizable shape ansatz for 
the equilateral shape that is often used in data analyses is the following [ST]: 

= "6 (J^ + 2 perm.) + 6 ^ + 5 perm.) - 12 , (6.22) 

and is shown in Fig. [5j As we can see, it represents the most important features of Fig. [3] & 
® 

The shape of S a is more complicated, but we expect they have the similar shapes as the 
equilateral one because their squeezed limits behave the same [50]. The three other shapes 
S e , S v and S s are all close to the local shapes as their squeezed limit scale as k\/ks for 
&3 < ki = k 2 . 

The scale dependence in P^, c s and A/S will introduce mild running for the three-point 
function. We usually regard only the contributions from c s and A/S as the running of the 
non- G aussianity. 

The underlying physics of the equilateral shape can be readily understood in terms of 
their generation mechanism. In single field inflation, the long wavelength mode that exits the 
horizon are frozen and can have little interaction with modes within the horizon. The large 
interaction only occurs among modes that are crossing the horizon at about the same time. 
These modes then have similar wavelengths. This is why the shape of the non-Gaussianity 
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peaks at the equilateral limit in momentum space. 

This physical origin also suggests the caveat that, as long as there are large interactions 
involving modes with similar wavelengths, an equilateral-like shape may arise. For example, 
such cases can happen in multifield models where there are particle creation [82j[83]. (See 
however [84|). 

• An example: Dirac-Born-Infeld (DBI) inflation. An explicit example of the 
above general results is the DBI inflation [5 ^ [71 |[85T - [92] . These inflation models describe a 
3+1 dimensional brane moving in warped extra dimensions. The location of the brane is a 
scalar field in 4d effective field theory, and it is the inflaton. The warped extra dimensions 
provide a non-trivial internal field space for the inflaton. In terms of the 4d effective field 
theory, the action is 



- / (fx. 



-fj 



f(<j>)- 1 y/l + Mg^d^rf- /(0)" 1 + V{4>) 



(6.23) 



The non-trivial part is the kinetic term involving the square- root. It can be understood as 
a generalization of the following two familiar situations. It is a higher dimensional general- 
ization of the action of a relativistic point particle 



(6.24) 



where the speed of light / x l 2 varies with x. It is also a relativistic generalization of the 
usual canonical kinetic term in the non-relativistic limit \f{(f))g^ u d ll (j)d v (j)\ <C 1, 



-g 



g^d^d u <P + V{<P) 



(6.25) 



Because the speed-limit of the inflaton f" 1 ^ 2 can vary in the internal space, if it can be made 
small enough near the top of potential where the inflaton is about to roll down, the warped 
space restricts the rolling velocity even if the potential is too steep for slow-roll inflation 
to happen. So the inflaton rolls ultra-relativistically, but with very small velocity, and this 
generates the DBI inflation. 

The physical consequence is now easy to obtain using the general results in this subsec- 
tion. In our notation the Lagrangian is 



P = -fVl-2/I + r l ~V . (6.26) 

The sound speed is 

c s = y/l - 2fX , (6.27) 
which is the inverse of the Lorentz boost factor 7, so c s <C 1. The component (16.18p vanishes 
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identically, and we have a large bispectrum of shape A c with size (16.191) . 

DBI inflation is still driven by the potential energy. The general single field inflation 
models also include the k-inflation [93], where the inflation is driven by the inflaton kinetic 
energy. Model construction of single field k-inflation can be found in Ref. [93~H95] . The 
bispectra for such models are computed in Ref. [50|l94] . 

Multifield generalization have been studied in Ref. [9614102] . where this type of kinetic 
terms are generalized to multiple fields. The three-point functions involving these different 
fields have the same or similar shapes. 

The current CMB constraint on the equilateral ansatz ( 16.221) is —214 < f c ^ L < 266 pp. 

6.2 Sinusoidal running: sharp feature 

Although various slow- variation parameters in ( 16 .5D have to be small most of the time during 
inflation, they can become temporarily large. Such cases can happen if there are sharp 
features in inflaton potentials or internal field space, so the behavior of inflatons temporarily 
deviates from the attractor solution, and then relaxes back within several Hubble time, or 
stay longer but with small deviation amplitudes. Motivations for such models include the 
following. It may be possible explanations for features in power spectrum [106-109], and if 
so the associated non-Gaussianity is a cross-check. And there are brane inflation models that 
are very sensitive to sharp features present in the potential or in the internal space [110] . 

As an example, we study a sharp feature in the slow-roll potential. The fact that a sharp 
feature in potential can enhance non-Gaussianities has long been anticipated and qualitative 
estimates have been made by different methods |103H105] . The precise method of analyzing 
the size, running and shape of such non-Gaussianities [SUES] is made possible with the 
developments of the formalisms that we reviewed in Sec. I3.1[ H] and 16.11 This will be the 
subject of this subsection. 

We start by studying the behavior of the slow-roll parameters. We use a small step in 
potential as an example and will ignore numerical coefficients. We use c ~ AV/V to denote 
the relative height of the step, and d the width of the step. In the attractor solution, the 
inflaton velocity is given by ~ V'/H ~ y/eV. As it falls down the step, the potential 
energy cV gets converted to the kinetic energy, so we have 

<P < s/V{c + e) . (6.28) 

The amplitude of density perturbations is given by P^ ~ H 4 /<j) 2 , so such a sharp feature 
causes glitches in the power spectrum. It will leave a dip with relative size AP^/P^ ~ 
\Jl + c/e — 1 since (f) increases first, followed by oscillations caused by a non-attractor com- 
ponent of the mode function before it settles down again in the attractor solution. To fit 
the CMB data, <fi cannot change much. As we can see, the sensitivity of the power spectrum 
to the step size c is proportional to e, and we need c/e < 1. Reducing the width d of sharp 
feature increases the amplitude of the glitches, but this is only for a large d over which the 
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inflaton spends more than one e-fold to cross. Further reducing d will not change the am- 
plitude of the glitches since (16. 28ft is saturated; but the sharpness will determine how deep 
within the horizon the modes are affected. 

So e does not change much, Ae ~ A((/) 2 )/H 2 ~ c. But it changes within a very short 
period, At ~ A0/0 ~ dj \/V{c + e). So 77 can be very large, 

Arj~ -£r ~ . (6.29) 

' HeAt de v ' 

It is also clear that the feature is associated with a characteristic physical scale and generates 
a scale-dependent power spectrum and higher order correlation functions. 

With these qualitative behavior in mind, we now study the three-point function. An 
important fact of the formalisms in Sec. H] and [6] is that the expansion is exact in terms of 
the slow-variation parameters. So it is valid even if these parameters are not always small, 
as long as the expansion in £ ~ O(10~ 5 ) is perturbative. 

In all terms in the cubic expansion (14. 8p . ( appears at most with one time derivative; the 
field redefinition gives a term that is proportional to rj at the end of the inflation; and the 
other terms are all suppressed by powers of e, which remains small even in the presence of 
a sharp feature. So the most important term is 

dtd 3 x^a 3 er]( 2 ( , (6.30) 

in which the coupling is proportional to 77. The correlation function (C(k 1 )C(k 2 )^(k3)) is 
dominated by 

\\uMend) j J ^ draper/ (ul(r)u* 2 (r)-^-u* 3 (r) + sym^) (2vr) 3 <5 3 (^ k*) + c.c. . (6.31) 

Precise evaluation of this expression has to be done numerically. But it is not difficult to see 
the generic properties of bispectra associated with a sharp feature. 

For long wavelength modes that already crossed the horizon at the time of the sharp fea- 
ture, kiT <C 1, the mode function is already frozen and the integration (I6.3ip gives vanishing 
contribution. For short wavelength modes that are still well within the horizon, the modes 
are not affected if their momenta are larger than the inverse of the time scale characterizing 
the sharpness of changes in slow-roll parameters. The modes most affected are those which 
are near the horizon crossing. These modes are all oscillatory, ~ e~ tkiT . As we have studied, 
7]' is temporarily boosted, so it can be roughly approximated as several hat-functions that 
satisfy j drrj' = 0. Examples of such behavior are shown in Fig. |6j If we simply approximate 
the hat-functions by several delta-functions, 7/ ex 8(r — r*), the integration (I6.3ip will give 



38 




Figure 6: Behavior of the slow-roll parameters for a step (solid line, c = 0.0018, d = 0.022Mp) 
and bump (dashed line, c = 0.0005, d = O.OIMp) sharp feature on |m 2 2 potential. Note 
that in the absence of the sharp feature, erf is of order O(10~ 4 ). 



something like 

S~/fefsin(^ + 0o) j (6 .32) 
where = 1 /t* is the scale corresponding to the location of feature, 0o is a phase and 

fig ~A V ~0 (^^) • (6.33) 

Comparing with the effect on the power spectrum, one can keep the size of glitches in the 
power spectrum small while make J'nl large, for example, by fixing c/e and decreasing d. 

This ansatz describes the most important running behavior of this bispectrum. Notice 
that the oscillatory frequency in the fc-space is of order 1/k*, which is the scale of the feature. 
A rescale in k* can be compensated by a rescale in all fcj. Also notice that the oscillatory 
frequency, 3/A;*, along the k\ = k 2 = A; 3 = k direction is 3/2 of that in the power spectrurd_|, 
2/K. 

1 For power spectrum, the sharp feature introduces a small non-Bunch-Davies component for the mode 
function. The oscillatory frequency in the power spectrum is determined by the phase of the coefficient for 
this component. This is obtained through matching conditions across the feature, and the phase is ~ 2fc/fc*. 
See Sec. 5. 3 of Ref. |110) for an example. 

From this result we can see that, observationally, while sharp features located at large scales (such as 
I ~ 30) introduce glitches that need to be distinguished from statistical fluctuations, those located at much 
shorter scales (such as I ~ 1000) introduce oscillatory modulation that coherently shifts all points over 
several acoustic peaks in the same direction, which is completely different from statistical fluctuations of 
data points. 
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Figure 7: Numerical result (dashed line) for the bispectrum running for a sharp step (c = 
0.0018, d = 0.022) along the k% = k 2 = k 3 = k direction, compared with the simple ansatz 
(solid line). 

In practice, ( 16 .32 j) is a crude ansatz that needs to be refined. First of all, we have only 
considered the modes that have not exited the horizon. For those did, as we mentioned, 
their correlation function is as small as usual. The ansatz needs to be cut off for the long 
wavelength modes K/k* <C 1. A more detailed analysis [BS] reveals, using the hat functions 
as an approximation of the slow-roll parameter behavior, that the bispectrum falls off as 
K 2 for these long wavelength modes. Secondly, the fact that in f 1 6 . 3 2 j) all short wavelength 
modes are equally affected is due to the sharp-change approximation. Smoother functions 
will only affect a finite range of modes within the horizon. So the amplitude of the ansatz 
should decay and how fast depends on the sharpness of feature. To take into account both 
effects, empirically, we can multiply ( 16 .32 j) with an envelop function 



where n and m are parameters chosen to fit the numerical results. For example, n = 2, 
m = 5 for Fig. [71 Lastly, in the very squeezed limit, k 3 <C ki\K/k#\, S can no longer be 
approximated as a function of K only and starts to have a non-trivial shape [65]. Here we 
concentrate on the signature running behavior. 
A numerical result with 



is shown in Fig. [71 A subtlety encountered in the numerical integration is how to handle 



oc {K/h) n e~ K/mk * , 



(6.34) 




(6.35) 
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the oscillatory behavior at r — > — oo. One can do a tilt into the imaginary plane, — r — > 
— oo(l + ie), as prescribed in the analytical procedure in Sec. 13.21 or more efficiently, perform 
integration by part to increase the convergence of the integrand at the r — > — 00 end. One 
may also try the method of Wick rotation, but this will first require solving the background 
equations of motion in the Wick-rotated space, since here we do not have the analytical 
expression for the mode function. 

Sharp features can also appear elsewhere instead of potentials, for example, in the internal 
warped space for DBI inflation |110j . The qualitative running behavior in bispectrum is 
similar, and overall large non-Gaussianities become a superposition of the approximate scale- 
invariant equilateral shape and the sinusoidal running. 

Non-attractor initial conditions can be included of sharp features, except that 

we only observe the relaxation part. 

6.3 Resonant running: periodic features 

In this subsection, we consider a different type of features. These features may or may not 
be sharp, but the most important property is their periodicity. Such features will induce 
an oscillatory component to the background evolution, in particular, to the couplings in the 
interaction terms. We denote this oscillatory frequency as uj. We know that each mode 
oscillates when it starts the life well within the horizon. This frequency keeps on decreasing 
as the mode gets stretched by the inflation, until it reaches H when the mode becomes 
frozen. So the mode scans through all frequencies that is larger than H, up to some very 
high cutoff scale such as Mp. Therefore, as long as 



the oscillatory frequency of the modes in the integral will hit uj at some point during the 
inflation. This cause a resonance between the couplings and modes, hence a constructive 
contribution to the correlation function [65J. Without the resonance, as we encountered 
previously, the highly oscillatory modes simply average out within the horizon. In contrast 
to the previous mechanisms, here the non-Gaussianities are generated when modes are sub- 
horizon. 

We now study the properties of such a non-Gaussianity, following Ref. |65j. 
To estimate the integral, we use the unperturbed mode function. Similar to the sharp 
feature case, we get 



uj > H, 



(6.36) 
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+ two perm. + c.c. . 



(6.37) 



41 



In this case, we are interested in the region \Kt\ 3> 1 in order to have resonance. So the 
last term dominates as long as the momentum triangle is not too squeezed so one of the k^s 
becomes < 1/r at the resonance point. The oscillatory coupling is dominantly contributed 
by rf. The integral is proportional to 

J drr sin (ut) exj>(iKr) . (6.38) 

This integral can be done analytically using the relation t « — H^ 1 ln(— Ht). But its most 
important properties can be understood as follows in terms of the physical picture that we 
described. 

First, let us look at its oscillatory running in i^-space. The phase of the background 
repeats itself after AN e = 2irH/u e-fold, during which the wave-number K changes by 
—KAN e . So the running of the non-Gaussianity in i^-space is also oscillatory with the 
period given by 

AK = KAN e = 2-kKH/uj . (6.39) 

Note that this period is changing with K in a specific way that we will see more clearly in 
a moment. 

Next, let us look at the size of the non-Gaussianity. Each i^-mode briefly resonates with 
the oscillatory coupling when its frequency sweeps through the resonance frequency u. Once 
its frequency differs from u by Au, the integration in the 3pt starts to cancel if is performed 
over At i ~ tt/Au. In the meanwhile it takes At2 ~ Au/(uH) to stretch the mode and 
change its frequency from wtow- Au. Equating At\ and At2 gives the time period over 
which the resonance occurs for this mode, 



(6.40) 



This corresponds to the number of oscillation periods 



uAt oj , 

that we need to integral over to estimate the resonance contribution. Note that one period 
in the integral (I6.38|) for K = u contributes tct*/K, where is evaluated at the resonant 
point. Multiplying the total number of the resonant periods ( I6.4ip . using the definition ( 15. 91) 
and Pq ~ H 2 /(87r 2 e), we see that the amplitude of S(ki, k 2 , k 3 ) is 

tsi ~ I? ~ ■ (6 - 42) 
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Figure 8: The running and shape of the resonance bispectrum ( 16.4511 with C = 100. 

Slow-roll parameters acquire small oscillatory components, and here rjA denotes the ampli- 
tude of such an oscillation. Other pre-factors of ki are cancelled according the definition of 
S and the S turns out to be a function of K only. In the last step of (I6.42p . we have listed 
the accurate numerical number, which differs from the estimate by a factor of y/2. 

Summarizing both the running behavior and the amplitude, we get the following ansatz 
for the bispectrum 

= InI sm (C HK/K)) , (6.43) 

where 

C = 2nK/AK = u/H (6.44) 

and fc* gives a phase. The argument Cln.fr in (16.431) appears because of (16.391) . This gives 
a scale dependent oscillatory frequency in the K-space. In fact, this kind of dependence 
makes the density perturbations in the resonance model semi-scale- invariant. We call it 
periodic-scale-invariant - they are invariant under a discrete subgroup of rescaling. Namely, 
the ansatz (I6.43j) is invariant if we rescale all fc, by nAK/ K = 2irnH/ u e-fold, where n is an 
integer. Other rescaling causes a phase shift. This property is a direct consequence of the 
symmetry of the Lagrangian. It is periodic, so invariant under a discrete shift of the inflaton 
field. This periodic-scale-invariance should also be respected by the full bispectrum results, 
as well as other correlation functions. 

As mentioned, we have derived this ansatz from the last term in (16.37)) . Other terms will 
become important in the squeezed limit. The full integration (16.37)) has been worked out 
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Figure 9: Numerical result (solid line) of the bispectrum running for the example (I6.46P 
(m = 3 x 10- 6 M P , c = 5 x 1(T 7 , A = 0.0007M P , « 15M P ), compared with the ansatz (RT33]1 
(dashed line). 



in [111] , and the leading order results are 

sin (C\n(K/K)) + i ^ ^ cos (Cln^/A;*)) + O ^ 



r res 

Jjvl 



(6.45) 



where 0(1/ C 2 ) terms are neglected because we need 1/C — H/uj <C 1 for large resonance. 
The numerical coefficient in ( I6.42p turns out to be y/ir/ (8a/2). As we can see, the extra terms 
satisfy the symmetry we mentioned and indeed give large corrections in the very squeezed 
limit, e.g. k 3 < kiH/u. These terms also ensure a consistency condition that we will study 
in Sec. 19.21 An example is plotted in Fig. [HJ The spike at the very squeezed limit is due 
to the second term in (16.451) . Overall, we see that the leading shape of this bispectrum is 
quite trivial, being almost a function of K only, until it gets to the very squeezed limit. The 
most distinctive feature of this type of non-Gaussianities is the running behavior captured 
in (I6.43p . Unfortunately, this ansatz is not factorizable if the .ff -range is too large. 

More arbitrary scale-dependence can be introduced if the features are applied over a 
finite range, or with varying periodicity and amplitude. 

As a useful comparison, the resonant running here and sinusoidal running that we stud- 
ied in the last subsection are clearly distinguishable from each other observationally. The 
resonant running oscillates with periods that are always much smaller than the local scale, 
AK <C K; the frequency has a specific scale-dependence, AK/K = const.; and the fre- 
quency in the power spectrum (~ sin(C In k) in /c-space) is exactly the same as that in the 
bispectrum (~ sin(Clnii') in K-space). In contrast, the bispectrum of the sinusoidal run- 
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ning oscillates with a fixed period that approximately equals to the scale at the location of 
the sharp feature, AK ~ the frequency is scale-independent; and the power spectrum 
(~ sin (2 A; /A;*) in fc-space) has twice an oscillatory frequency of the bispectrum (~ sm(K/k*) 
in .fT-space). 

As an illustration, we look at an example, 



1 



V(cf)) = -m 2 2 



1 + c sin — 
A 



(6.46) 



In this example, the inflaton is rolling over the small but periodic ripple laid on the potential. 
This induces an oscillatory component in the slow-roll parameters with an amplitude tja ~ 

\^Qcm(J)/A 2 and a frequency u ~ 0/A 2m/(v / 6A). So we have 



res 
NL 



A 5 /201/2 



(6.47) 



and 



C « 2/(0A) . (6.48) 

A numerical example is shown in Fig. E0 As we can see, the ansatz (I6.43P gives a very 
accurate fit to the actual running behavior. The mode function and power spectrum are the 
superposition of the usual unperturbed solution and a small oscillatory component [65 |I1121 
1113] . We can choose parameters so that the size of the ripples on the power spectrum is 
small, but bispectrum is made large. This is because the non-Gaussianities rise more quickly 
if we increase the frequency, while the mode function has difficulty responding efficiently 
when the external source oscillates too fast. This mechanism may be realized in terms 
of brane inflation [110J where the periodic feature comes from duality cascade in warped 
throat |114j . or the monodromy inflation [115(1116] where the periodic feature comes from 
instanton effects |112|lll3j . 



6.4 Folded shape: a non-standard vacuum 

In this subsection, we study the effect of non-standard vacuum on the primordial non- 
Gaussianities. We consider a different wave-function from the Bunch-Davies vacuum when 
modes are well within the horizon. To start, let us first discuss several motivations for this 
case. 

• A non-Bunch-Davies vacuum can actually occur much more simply than it might sound 
like. Any deviation from the attractor solution of the inflaton generically generates 
a component of non-Bunch-Davies vacuum. This is because a general mode function 
is a superposition of two components, ci(k)u(k, t) + C2(k)u*(k, t), and in attractor 
solution we choose one of the component asymptotic to the Bunch-Davies vacuum. 
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A disturbance will generically introduce a mixture with the other component. In 
this sense, we have already encountered such a case when we studied the effect of 
a sharp feature in Sec. 16.21 Indeed, after the inflaton crosses the sharp feature, the 
oscillatory behavior in the power spectrum is precisely due to the superposition of the 
second non-Bunch-Davies component for some finite /c-range. For an infinitely sharp 
change, such a disturbance with a small amplitude extends to all k that have not exited 
the horizon at the time of sharp feature. An analytical illustration can be found in 
Sec. 5.3 of Ref. [110] . The location of the sharp feature can become superhorizon at 
the present time, but its influence has extended to much smaller scales and becomes 
observable. The resonance case in Sec. 16.31 is another type of example. An analytical 
illustration can be found in Sec. 3.3 of Ref. [112] . For non-Gaussianities studied in 
Sec. 16.21 & I6.3[ we only concentrated on the effects caused by the features in slow-roll 
parameters. In analytical analyses, we approximated the mode function by the Bunch- 
Davies component and ignored the disturbance. The study of this subsection can be 
regarded as the complementary analyses on the effect of a different mode component. 

• In inflationary background, modes can be quantized in terms of time- dependent cre- 
ation and annihilation operators, ak(t) and a_ k (i). The Bunch-Davies vacuum is de- 
fined as the vacuum annihilated by a^(t) as t — > — oo. If a different adiabatic vacuum is 
defined which is annihilated by cik(^o) & t a finite to, we introduce a non-Bunch-Davies 
component. For example, see |117[I118] . The origin and magnitude of such a com- 
ponent have been debated and studied by many papers, often under the name of the 
"trans-Planckian effect". See Ref. |119[I120] for summary and references. 

• There are inflation models where the scale of new physics can be very low. In particular, 
in warped space it is proportional to the exponentially small warp factor. In some DBI 
inflation models [86, 92], the speed limit of the inflaton and the scale of new physics are 
both related to the warp factor in such a way that the local warped new scale can drop 
near or even below the Hubble energy scale in certain epoch of inflation. Clearly the 
simple scalar field Bunch-Davies vacuum is no longer sufficient. Such models further 
open up the possibilities of vacuum choices. 

After these discussions, let us now focus on a specific simple problem [50]. We modify 
the wave-function of the Bunch-Davies vacuum by a small second component and examine 
its consequence for the three-point function calculated in Sec. 16.11 We consider the general 
single field inflation with a small sound speed c s or a large A/E [50 |I121] . 

So the mode function is 



where |C_| <C 1 and can be /c-dependent. In the first example above, the extra component 
starts at a specific time in the past. In the second class of examples, it may start either at 




iH 



(6.49) 
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Figure 10: Shape of S\ (truncated). 



a specific time or specific energy scale. To see a common feature without addressing these 
model-dependent issues, we look at the simple limit where the r in (I6.49P can go all the way 
to — oo. The computation of the correlation function is essentially the same as in Sec. 16.11 
The leading order correction to the bispectra is obtained by replacing one of the three Mfc(r) 
in the integrand by its C_ component. So it simply replaces one of the k^s in the shapes 
with —hi. For example, the correction to S\ is 

5 / 1 2A\ 3hk 2 k 3 

A = ~ S ) — 2 — 

x f ^ML + + , gg(%M ) . (6 . 50) 

The shape of S\ is shown in Fig. [TUJ The most important feature of this shape is 
the enhancement at the folded triangle limit, e.g. k\ + k 2 — k 3 = 0. The detailed form 
of enhancement is model dependent. For example, it is different for another shape S c . 
The divergence in this folded limit occurs due to our simple limit of taking r to — oo. 
Imposing some kind of cutoff at the lower limit of r will eliminate this divergence, although 
as mentioned the detailed modification will be highly model dependent. For example, a 
simple constant cutoff r c will introduce a factor of 1 + (\x 2 c — 1) coss c — x c sina; c for each 
of the three terms in ( 16.501) . where x c = (k\ + k 2 — kz)c s r c or its cyclic. Very close to the 
folded limit, i^c s |r c | <C 1 (K = k\ + k 2 — k 3 or its cyclic), this regulates away the divergence; 
away from the folded limit, ifc s |r c | 3> 1, these extra factors are unity on average but with 
oscillations. These oscillation can be either physical, or regarded to be zero if x c is within 
a regulation scale which exists since the non-Bunch-Davies component is present for a finite 
time in the past. 
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(a) (b) 

Figure 11: Two ansatz for the folded shape, (a) Eq. (16.511) . (b) Eq. f !6.52j) with k c /k\ = 0.1. 

The case for slow-roll inflation is qualitatively similar, and more examples of the bis- 
pectra shapes and the observational prospects are discussed in |122[I123] . In this case, the 
proportional parameter for the bispectra amplitude is no longer enhanced by l/c% or A/S, 
but < 0(1). 

In order to facilitate the data analyses, a simple ansatz has been proposed in Ref. [123J, 
S^Ui = 6 (J^ + 2 perm.) -6^ + 5 perm.) + 18 , (6.51) 

which represents certain important features of this kind of bispectra. It has a smooth rising 
behavior in the folded limit. This ansatz is plotted in Fig. UTT a). Since the real shape has a 
model dependent cutoff, it remains open questions how sensitive this is to data analyses and 
how well the ansatz (I6.5ip represents it. We can also write down an ansatz which is more 
directly motivated from the example (I6.50p and the comments after that equation, 

= + 2 perm ' ' (6 - 52) 

where the cutoff scale k c = l/(c s r c ) is a parameter. For k% + k 2 — k 3 ^> k c and cyclic, we 
have neglected the oscillatory part and only taken the average. In this ansatz we can change 
the powers in the numerator and denominator to model model-dependent variations. The 
relation (k c + h + k 2 - k 3 )~ n = (r(n))" 1 / °° dt t n - 1 e ^ kc+kl+k2 - k ^ t may be used to factorize 
the ansatz. This ansatz is plotted in Fig. ITTT b). 

Another type of non-Bunch-Davies vacuum, namely an n-particle state built on the 
normal Bunch-Davies vacuum, was studied in |124[ll25j and the non-Gaussianities were 
found to be unobservable. 
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7 Quasi-single field inflation 



Having considered single field inflation, we now relax the condition on the number of fields. 
At least during inflation, we only need to consider quantum fluctuations of light fields, since 
if the mass of fields are very heavy, (here the relevant scale is m > H), they contribute 
only classically and determine the classical inflaton trajectory. Multiple light fields can 
arise naturally if we consider the inflation models as the consequence of a UV completed 
framework. However, as discussed in Sec. I5.1[ due to the back-reaction from the inflationary 
background, the mass of light fields are naturally of order H. The potential with such a 
shape is too steep for slow-roll inflation. 

Therefore, as a natural step beyond the single field, let us consider slow-roll models with 
one inflationary direction, and one or more other directions that have mass neither much 
heavier nor much lighter than H. We will call the quanta in the inflationary direction as 
the inflaton and its mode the curvature mode, and the others isocurvaton and isocurvature 
modes. We call these models the quasi-single field inflation models [126, 46J. 

Note that the thematic order in this review is not chronological. The non-Gaussianities 
in this type of models were not computed until very recently for a couple of reasons. If the 
mass of particles is of order O(H) or larger, the amplitude of these fields decay exponentially 
in time after horizon-exit. So they would not seem to be important for super- horizon pertur- 
bations even if they couple to the curvature mode. As we will see, however, their amplitudes 
at or near the horizon-exit are enough to make them interesting. What really suppresses 
their contribution is the fast oscillation behavior present for m 3> H . Methodologically, 
isocurvature-to-curvature transition for non-Gaussianities was studied restrictively in the 
regime of super- horizon classical evolution in multi-field space |127H134] . which we shall 
explain in more details in the next section. However, for quasi-single field inflation mod- 
els, a full quantum computation in the in-in formalism is necessary to properly include the 
contributions from both the horizon exit and the superhorizon evolution. 

7.1 Intermediate shapes: massive isocurvatons 

There are potentially different ways massive isocurvatons can be coupled to the inflaton. We 
currently do not have a general approach in terms of model building. So what we shall do 
is to first study this problem through a simple example, and then discuss the features of the 
results that can be regarded as generic signatures of this class of models |126[ |4"6]. 

We consider the case where the inflaton is turning constantly by going around (a fraction 
of) a circle with radius R in the angular 9 direction. See Fig. [T2J All the parameters, such 
as R and couplings, are assumed to be constant during the turning. We call this assumption 
the constant turn case. In the 9 direction the potential is the usual slow-roll potential V sv (9). 
The field in the radial direction is denoted as o and has mass of order H, and lifted by the 
potential V(o~). For such a turning trajectory, it is convenient to write the action in terms 
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Figure 12: Quasi-single field inflation with turning trajectory. The field 6 and a are in the 
polar coordinates. The 9 is the inflationary direction with a slow-roll potential. The a is the 
isocurvature direction, which typically has mass of order H. 



of fields in the polar coordinates, 6 and a, instead of in the Cartesian coordinates, 



Sn 



d A X\J—g 



-~{R + afg^d^dj - l -g^d,ad u a - V BT (6) - V(a) 



(7.1) 



The potential V{a) balances off the centrifugal force necessary for the turning and traps the 
field at the bottom of the effective potential, V e s(a) = — \Qq(R + cr) 2 + V(a). We define the 
minimum of this effective potential to be a = 0. We expand the effective potential as 



V c fi = const. -| — [V 



(7.2) 



where 9 is the turning angular velocity and the primes on V denote derivatives with respec- 
tive to cr. 

To study the perturbation theory, we perturb the fields in the spatially flat gauge, 



0(X, t) = Q (t) + <J0(X, t) , a(x, t) = 5cr(x, t) , 

and obtain the following Hamiltonian, 

H = a 3 \\R 2 5e? + |(^,) 2 + isa/ + ^(dMi) 2 + \m 2 5a] 



(7.3) 



(7.4) 



= —c 2 a 3 5ai59j , 
U\ = c 3 a 3 5a 3 j , 



(7.5) 
(7.6) 
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where 

c 2 = 2R9 , c 3 =lV", m 2 = V" + 76 2 Q (7.7) 
o 

are all constants. Terms suppressed by (9(e) have been ignored in this gauge. The curvature 
perturbation ( is most transparent in another gauge, the uniform inflaton gauge, where 

0(x, t) = 9 (t) , <r(x, t) = (To (t) + 5a(x, t) , (7.8) 

and the spatial metric is 

Mx,*) = aW C(x '% • (7-9) 

In this gauge, ( appears in the metric as the space-dependent rescale factor and the fluctu- 
ations in the inflaton is shifted away. The relation between ( and 69 is the gauge transfor- 
mation. At the leading order this is 

C « -5-60 . (7.10) 

We will calculate the correlation functions in terms of 59 and then use this relation to convert 
them to those of (. The full perturbation theory that one obtains in the uniform inflaton 
gauge justifies the above omission of several 0(e) terms in the spatially flat gauge [46J. 
There are several important points for this Hamiltonian. 

First, the kinematic Hamiltonian (17.41) describes two free fields in the inflationary back- 
ground. One is massless and has the familiar mode function, 

u k = H — (l + ikT)e- tkT . (7.11) 
RV2k* 

Another is massive and the mode function is 



vy- = —le 

2 



3 i("H)5^#(- r )3/2#(i)(-jfe r ) , (7.12) 
where 

v = y/9/4- m 2 /H 2 . (7.13) 

For < m/H < 3/2, the amplitude of the mode decays as (— r)~ u+3 / 2 after horizon-exit 
kr — > 0. The lighter the isocurvaton is, the slower it decays. At the v — > 3/2 (i.e. m/H — > 0) 
limit, the amplitude is frozen. For m/H > 3/2, v becomes imaginary, the mode not only 
contains a decay factor (— r) 3//2 but also an oscillation factor r v . This oscillation is marginal 
for m ~ H, but if m ^> H, it causes cancellation in the integrals of the correlation function 
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Figure 13: Feynman diagrams for the transfer vertex (a), corrections to the power spectrum 
from isocurvature modes (b), and the leading bispectrum (c). 

and is equivalent to factors of Boltzmann-like suppression ~ e ~ m / H \ We will consider the 
case < v < 3/2. 

Second, there is a sharp contrast between the V™ for the slow-roll inflaton field and 
the V" for the massive field a in the non- inflationary direction. The former has to be 
very small, ~ 0(e 2 )P^ 2 H , in order to maintain the smallness of the slow-roll parameters. 
(Here we use e to denote collectively all slow-roll parameters, e = —if /if 2 , rj = e/eH, 
and £ = r]/r)H.) Consequently it contributes 0(e 2 ) to the /nl of bispectrum in slow-roll 
inflation, generally smaller than the 0(e) contributions from the other terms in the same 
model. However, for quasi-single field inflation, the direction orthogonal to slow-roll does 
not have to satisfy the slow-roll conditions, and V" is almost unconstrained. For example, 
in the inflationary background, it can be of order H; and similarly, V"" can be of order one, 
etc. This isocurvaton self-interaction f lT.6 j) becomes the source of large non-Gaussianities. 

Third, the coupling between the isocurvaton and inflaton appears as a form of a two-point 
vertex operator in ([7.5]) . We treat this term as part of the interaction Hamiltonian, and it is 
represented by the transfer vertex in Fig. [12] (a). The strength of the coupling is determined 
by the turning angular velocity 8q in this model. This coupling is responsible for the trans- 
formation of the isocurvature perturbations, in particular their large non-Gaussianities, to 
the curvature perturbation. 

We calculate correlation functions corresponding to the Feynman diagrams Fig. [T3] in 
terms of the in-in formalism, which we reviewed in Sec. 13.11 As an illuminating example to 
illustrate the different advantages of the three forms of the in-in formalism, we recall from 
Sec. 13.11 that the three-point function can be written in the following forms. The original 
definition (I3.44[) and (I3.40p . which we refer to as the factorized form, leads to 



12^3^(0)^(0)^3(0) 



■n 



x Re 



/ df x a 3 (fi)u^(7i)i4(fi) 



dr 2 a A (T2)v pi (f 2 )v P2 (T2)v K (T2) 



— oo 
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X 



PT\ 

dn a^r^v* (t^u* (n) / dr 2 a 3 (r 2 )v* Jr 2 )u'* (r 2 ) 



x (27r) 3 5 3 (^^ pj) + 9 other similar terms 

i 

+ 5 permutations of p» . 



(7.14) 



The perturbation theory here starts from the fourth order. The reorganized commutator 
form (I3.4ip leads to 

(59 3 ) = 12c 3 C3U Pl (0)u p2 (0) Mp3 (0) 



x Re 



dn / dr 2 / dr 3 / dr 4 ]"J (a 3 (r;)) 

-oo J —oo J —oo J —oo 

x a(r 2 ) (u' Px {tx) - c.c.) (n)v^ (r 2 ) - c.c.) 

x (%^T 2 )^ 3 (T 4 )Up* (r 4 ) - c.c.) u P2 (r 2 )u* 2 (r 3 )«p 2 (r3 

x (27r) 3 (5 3 (^^ + 2 other similar terms 

i 

+ 5 permutations of p» . 



(7.15) 



In the IR (r — >■ 0), each of the ten terms in the factorized form diverge as r 3 &u for 
3/2 > v > 1/2 (0 < m < \/2H); while in the commutator form, various subtractions off 
the complex conjugates and the requirement that the final result has to be real make such 
divergence explicitly disappear. 

In the UV (r — > — oo), each factor of the multiple integral that integrates from — oo to 
has a definite convergent direction if we choose one of the two contour tilts, T{ — > — oo(l±ie), 
accordingly. Or more efficiently, by a Wick rotation — > ±izi. This would have been the case 
for the commutator form if we can break up the integrand into individual terms. However 
in order to achieve the explicit IR convergence, as we saw above, these terms have to be 
grouped; but then they have contradicting convergence directions. 

To take advantage of both forms, we introduce a cutoff r c , and write the IR part (r c < 
t < 0) of the integrals in terms of the commutator form, and the UV part (r < r c ) in terms 
of the factorized form, in the following mixed form (13 .42 p . 



dTi {commutator form} / dr i+ i 



dr n {factorized form} . 

(7.16) 

This shows explicitly both convergence behavior of the correlation function. Combining with 
Wick-rotations of the integration contours in the UV, this form provides an efficient way to 
evaluate the correlation functions numerically. The shapes of bispectra are presented in 
Fig. M for v = 0,0.3,0.5,1. 
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Figure 14: Numerical results for the shapes of bispectra with intermediate forms. We plot 
S with v = 0, 0.3, 0.5, 1. The plot is normalized such that S = 1 for p 1 = p 2 = p^- 

To better understand the shapes analytically, we can work out the squeezed limit (p^ <C 
Pi = P2) of the three-point function, 

<M 3 > ~> ft 7 u l a+v s(u) (2nmJ2 Pi) > ( 7 - 17 ) 
Pi P2PI i 

where s{u) is a ^-dependent numerical number. 

Recall that the squeezed limit of S for the equilateral shape goes as pz/pi, while for the 
local shape (ps/pi) -1 - Both the numerical results in Fig. [T4] and the analytical results in 
Eq. ( I7.17P show that here we have a one-parameter family of shapes, ~ {j>zlPi) x ^ 2 ~ l ', he 
between the two. We call them the intermediate shapes. 

The physical origin of such shapes can be understood as follows, and should be a generic 
signature for the quasi-single field inflation models. As we have seen, the large equilateral 
non-Gaussianity arises because the interacting modes cross the horizon around the same 
time. The shape of bispectrum peaks at the equilateral limit where the modes all have 
comparable wavelengths. As we will see in Sec. El the large local non-Gaussianity arises due 
to the classical non-linear evolution of superhorizon modes in the multifield space; so the 
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Figure 15: Shape ansatz (17.181) for the intermediate forms. 



interactions are causally disconnected and behave local in position space. This is non-local 
from the momentum space point of view. So the shape of bispectrum peaks at the squeezed 
limit. Now for quasi-single field inflation, the large non-Gaussianities come from the massive 
isocurvaton. Depending on the mass, these modes either decay right away after they exit 
the horizon (for m > y2H), or survive for a long time at the super-horizon scales (for 
m < \/2H). In the former case, the generation and transfer of non-Gaussianities maximize 
for modes that are exiting the horizon around the same time, resulting in quasi-equilateral 
shapes; in the latter case, the generation and transfer of non-Gaussianities happen in a 
superhorizon fashion, resulting in quasi-local shapes. In this regard, let us look more closely 
at the special limit m/H — > [y — > 3/2). 

In this massless limit, an infrared cutoff to the integrals are necessary. Otherwise the 
transfer will last for ever for the constant turn case. The cutoff corresponds to the ending 
of the turning. Let us discuss the following two cases. First, we still keep V" large. Our 
analyses still apply in this case. Interestingly, the shape of the bispectrum goes to that 
of the local form in this limit. As we will explain in Sec. 18.11 this is a generic signature 
of a massless isocurvaton. The infrared e-fold cutoff will introduce some running in the 
/atl because different modes experience different turning e-folds. Second, we would like to 
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make the isocurvature directions flat so this becomes a two-field slow-roll inflation models. 
Such models were intensively studied and it is known that the isocurvature modes can be 
transferred to the curvature mode by turning. However, since V" ~ 0(e 3 ^ 2 )H 2 /Mp is 
required to maintain the small slow-roll parameters, the contribution we computed here 
generates too small non-Gaussianity. We expect contributions from other terms are small 
as well. So it is much more difficult to generate large non-Gaussianities in such models, 
essentially because imposing the slow-roll conditions in all directions are too restrictive. 

To connect with data analyses, guided by the numerical results and analytical squeezed 
limit, we can use the following ansatz to describe the full family of shapes, 



,int _ 3* fTM +P2+ P 2 3 )(PlP2P 3 )l- 




cunt _ " JMLKfL 1 ^2 1 fj)\fVfAfi)- /r, -I q\ 

° ansatz 1n , . 7 o ■ l'- ic V 

10 (Pi + P2 + Ps)^ 3V 

These shapes are shown in Fig. [T5j Comparing with Fig. [TU we can see that they match 
quite well except near v = 0, around which it is better represented by another ansatz in 
Ref. |46j. The size of this bispectrum is 



(7.19) 



where ~ 6.1 x 10 -9 and 8q is the turning angular velocity. The ol(v) is a positive numerical 
number which, depending on u, can give an additional enhancement factor of order Nf (Nf 
is the turning e-folds). Since (8q/H) 2 and V"'/H are the expansion parameters in the 
perturbation theory, they have to be small to trust our calculation. Nonetheless this is not 
the model-building requirement. 

The fluctuations of more massive (m > 0(H)) fields may become important if they play 
a role later in the reheating |135[I136] . Such cases typically require some tunings for special 
conditions, so that the highly suppressed fluctuation amplitude can become important. 



8 Multifield inflation 

As we have seen in Sec. I7.1[ if we take the isocurvaton mass to zero in quasi-single field 
inflation while keep the nonlinear self-couplings of the isocurvaton V" large, the shape of 
the large bispectrum in the squeezed limit approaches the local form. The local form is in fact 
the earliest and most well-studied example of non-Gaussianities |127[I137[I138[ I(36]. although 
it was first found to be small as we have seen in Sec. HI As we will explain in this section, a 
large local form is a signature of massless isocurvatons that have large non-linear evolution in 
multifield space. We have arrived this shape from the in-in formalism by taking the massless 
limit. But if we stay in this limit, there is an easier formalism, the 5N formalism |139H141] , 
in which the underlying physics of the local shape becomes transparent. 
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8.1 Local shape: massless isocurvatons 

We recall that, in single field inflation, if we use the uniform inflaton gauge where there are 
no fluctuations in the inflaton field, the scalar perturbation ( enters in the scale factor as 
a 2 e 2C - . For superhorizon modes, ( is frozen. If we look at the different comoving superhorizon 
patches, they are causally disconnected from each other. So they evolve independently and 
locally in space. In such a gauge, the only difference is a space-dependent scale factor. This 
is also called the separate universe picture. The primordial curvature perturbations manifest 
themselves as the different number of expansion e-fold, 5N, at different positions. 

We would like to generalize this picture to the multifield case in the following 8N formal- 
ism. We will resort to a simple version of 5N formalism stated below, which is of course a 
consequence of the in-in formalism, but formulated from a simple perspective which clearly 
illustrates the points in this section. Otherwise, as we will explain, in the most general sense 
the 5N formalism should be defined as the in-in formalism written in terms of specified 
gauges. 

• We consider a set of scalars fa during inflation. Inflaton is one of them but can be 
different linear combinations of 0j's as a function of time, and the other orthogonal 
fields are called the isocurvatons. All the modes that we are eventually interested in 
should all have become superhorizon when the initial slice (specified below) is chosen. 
We look at different horizon-size patches and label them with the coarse-grained co- 
moving coordinate x. In the in-in formalism, the superhorizon modes behave as the 
c-number time- dependent background for each comoving patch. So we evolve these 
patches independently and classically. 

• We pick an initial spatially flat slice, on which there is no scalar fluctuations in the 
metric and all the fluctuations are in the scalar fields 0o« + <5<M X )- We assume that we 
know the statistics of such fluctuations. 

• We pick the final uniform density slices. Relative to the unperturbed and perturbed 
initial spatially flat slices, we have, respectively, the unperturbed and perturbed final 
uniform density slices. For single field inflation, these two final surfaces are the same. 
For multifield models, they are generally different. Such final slices have the properties 
that the universe has the same energy densities and field configurations everywhere on 
them. They can be chosen during either the inflation or the reheating. After that, 
every separated universe will have the same evolution. The only difference is the scale 
factor. This is the analogy of the uniform inflaton gauge in single field inflation. We 
study the cases where such slices exist. 

• We evolve the unperturbed 0oi m the initial slice classically to the unperturbed final 
slice, and denote the number of e- folds as N ((f) 0i ). This is the unperturbed e-fold 
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number. We evolve the perturbed 0oi + 50;(x) m the initial slice classically to the per- 
turbed final slice, and denote the number of e-folds as iV(0 Oi + 5</>i(x)). The difference 
between them 



5N = N(<j )0l + 6Mx))-No(M , 



(8.1) 



is the curvature perturbation £. Here iVo is a constant that can be shifted to make 
(SN) zero. 

• We expand 



where the subscripts on iV denote the partial derivatives with respect to 0j. For ex- 
ample, Nj = (dN/d(f)i)(dN/d(f)j). Repeated indices are summed over. The correlation 
functions of ( can then be computed as the classical averages of the products, such as, 



• We have assumed that the statistics of the 50* (x) are known on the initial slice. But 
this is not always easy to get. So we will consider the simple cases where this statistics 
can be approximated as Gaussian. Otherwise, calculating such initial statistics requires 
using the full quantum mechanical in-in formalism. 

Most generally, one identifies SN with the scalar curvature ( in the uniform infla- 
ton gauge; and the relation between SN and 8<f>(x) in the SN formalism is the gauge 
transformation between the uniform inflaton gauge and the spatially flat gauge. Cal- 
culating the correlation functions for ( becomes calculating those for 50 (x) using the 
in-in formalism. An example is the one we have seen in Sec. 17.11 

• So far we have not used the condition that the isocurvatons are massless (m <C H). If 
they are massive, after horizon exit the modes decay. So in the superhorizon classical 
regime, where the SN formalism is supposed to be useful, we are back in the single 
field inflation. Sub and near horizon perturbations should be computed by the full 
in-in formalism. Therefore having massless isocurvatons opens up classical multifield 
space in which we can have sizable SN defined in ( 18. II) . 

Now let us consider the Gaussian fluctuations 5 fa. From Sec. [2, we know that for massless 






scalars 
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where if* is the Hubble parameter when the corresponding mode exits the horizon. If the 
scalars are not exactly massless, H will have a running dependence on k\ caused by the 
decay of the amplitude. Using 

5N(k) = iW,,(k) + -Ny J 7^i^(k - k')^(k') , (8.6) 

we get the power spectrum 

2kf 



(C(k 1 )C(k 2 )) = (2 7 r) 5 ^5 3 (k 1 + k 2 ) , (8.7) 



where 



and the bispectrum 



(C(k 1 )C(k 2 )C(k 3 )) = NqNttyS- + 2 perm.) (2 7 r) 3 5 3 (^k l 



.9) 



According to the definition (I5.9p . the shape function is 



~<\oc ° ^loc / ^1 



S -i5^lfefc + 2perm | ' (810) 



where 



fl oc_ 5 NyNjNj 
Tnl ~ 6 (Np 2 ' 



As usual, we have ignored a mild running from the power spectrum P^. The shape (18.101) is 
called the local shape, which we plot in Fig. [16j This form is already factorizable. 

The physics of this shape can be understood from the derivation above. As explicitly 
demonstrated in (18. lft - (18. 4p . this non-Gaussianity is generated locally in position space for 
superhorizon modes. After Fourier transform, it becomes non-local in momentum space. 
That is the reason that the shape peaks at the squeezed limit. 

If the perturbation <50(x, t) on the initial spatially flat slice cannot be approximated as 
Gaussian, the shapes of final bispectra can be more complicated. 

• An example: the curvaton model. We use the curvaton model |143H151] as an 
example to illustrate the generation of large local non-Gaussianity. We also use it to demon- 
strate the 6N formalism. 

In this model, we assume that during inflation there is another light field o with the 
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i.7 - 

Figure 16: Shape of the local form (18.101) . 

potential 

V{a) = , (8.12) 

and m a <C H. This field is called the curvaton field for reasons that will become clear 
shortly. The energy density of the curvaton field is negligible initially. During inflation, it 
fluctuates and obtains the primordial amplitude ex* = <to + 5<t(x), where * denotes its value at 
the horizon-exit and after that the amplitude is approximately frozen. These perturbations 
are Gaussian for the potential (I8.12p with the canonical kinetic term, but can be more 
complicated otherwise. Here we study the simple Gaussian case. After inflation, these a- 
modes remain frozen until the Hubble parameter drops below m a . Then the a- field starts 
to oscillate around the bottom of the potential and behavior as matter. The Universe is still 
radiation-dominated. The fraction of the matter energy density quickly grows, because the 
matter dilutes as a~ 3 while radiation a~ 4 . The a- field decays to radiation when it reaches 
its lifetime. 

Another assumption of the curvaton model is that the primordial fluctuations in the 
inflaton field is much smaller than what is needed to achieve ( ~ O(10~ 5 ), although their 
total energy density may still be the dominant one. So the primordial curvature perturbation 
is contributed by the fluctuations in the a field, hence the name curvaton field. 

At the initial spatially flat slice t , we denote the radiation and curvaton density as p r o, 
p m o, respectively, and the scale factor as ao- Both components initially redshift as radiation. 
This lasts until the Hubble parameter reaches m a at t 1 . We denote the scale factor at ti as 
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a\. The Friedman equation at t\ is 

ZMlml = (^j ( Pr0 + p m0 ) « Pro • (8.13) 

After this the curvaton starts to oscillate and behave as matter. Denote the decay rate of 
the curvaton as T. We use the sudden decay approximation and assume that they decay 
instantaneously at the epoch H = T, because a process that has the time scale T falls into 
the Hubble expansion epoch with H = 1/T. We denote this instant as t 2 . The Friedman 
equation at t 2 is 

3M 2 T 2 = M 4 pro + M 4 3 Pmo . (8.14) 

Because at t 2 the universe has the same Hubble parameter (hence the same energy density) 
everywhere, this is the final uniform density slice. After that, both components become 
radiation and the evolution everywhere is the same. So we want to work out the expansion 
efolds N from to to t 2 as a function of the initial field value a. From (18. 14ft . we get 

e _47V + e -3N a = congt ^ ( gl5 ) 

where a = (a /ai)(p m o/p r o)- From (18. 131) we can solve for ao/ai which is independent of o 
at the leading order, and from ( 18 . 1 2[) we know that p m Q is proportional to a 2 . Therefore a 
is proportional to a 2 . Also note that r p = (a 2 /a )a = e N a = p m / pr\t 2 is the ratio of the 
energy density between the curvaton and the rest of the radiation at t 2 . Using these simple 
facts, we can differentiate (18. 15ft with respect to a once and twice, and get 

loc _ 5 AW _ 5 5(4 + 9r p ) 
/iVL-g^T-^r 12(4 + 3r p ) • [ * Ab) 

In terms of the definition r = 3p m /(4p r + 3p m )| t2 = 3r p /(4 + 3r p ) often used in the literature, 

/m-h-l-T- (8 - 17) 

So Large local non-Gaussianity arises if r 1. Note that although ( 18. 15ft only depends 
on a, this is a multifield model because the curvaton takes effects during the reheating. In 
some simple models in which the curvaton leads to non-adiabatic perturbations between 
dark matter and photons, r is tightly constrained by observations |142j . 

The large local form has been studied most extensively in the past. Variety of possibilities 
exist. They all share the common feature that non-Gaussianities are generated by some 
massless isocurvaton fields which acquire the superhorizon evolution during the inflation. 
For example, in multifield slow-roll inflation a turning trajectory [152J can transfer non- 
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Gaussianities from other directions to the inflaton direction [127H134] . But it is found to 
be very difficult to make non-Gaussianities large essentially because the very restrictive 
slow-roll conditions in all directions. In modulated reheating |153[I154] or preheating |155f - 
159J scenarios, the role of isocurvatons are played by the massless fields which control the 
couplings during the reheating or preheating. Thus they create a large local non-Gaussianity 
in a similar fashion as the curvaton model. 

Local form is also found in different contexts, such as models with special types of mas- 
sive gauge fields that acquire super-horizon evolution [160H162| . some non-local theories of 
inflation |163[ll64j . and certain cyclic universe scenario as alternatives to inflation [165-169J. 

The current CMB constraint on the local bispectrum is —10 < /]^£ al < 74 pp. Current 
constraint from large scale structure gives —29 < /^£ al < 70 [170] . Variety of methods have 
been invented to measure the local and other different forms of non-Gaussianities [1711 - 1184] . 

9 Summary and discussions 
9.1 Summary 

In this subsection, we summarize the main results of Sec. [6] - El Non-Gaussianities, con- 
ceptually being the expectation values of perturbations in a time-dependent background, 
are defined by the first-principle in-in formalism. Physically, having large primordial non- 
Gaussianities means that there are large non-linear interactions of some sort determined by 
certain dynamics during inflation. Measuring them tells us the nature of the dynamics. 

• Equilateral shape and higher derivative kinetic terms. 

In single field inflation, the long wavelength modes that already exited the horizon 
are frozen. They cannot have large interactions with short wavelength modes that 
are still within the horizon. When modes are well within the horizon, they oscillate 
and the contributions to non-Gaussianities average out. Therefore the only chance to 
have large interaction is when all modes have similar wavelengths and exit the horizon 
at about the same time. Theories with higher derivative kinetic terms provide such 
interaction terms. This is why the resulting bispectrum shape peaks at the equilateral 
limit in the momentum space. It drops to zero at the squeezed limit k 3 <C h\ = ki as 
ks/k\. It happens that, when these higher derivative terms become important enough 
so that the inflationary mechanism is no longer slow-roll, these large non-Gaussianities 
become observable. The forms of the bispectra are given in Eq. (16.161) and (16.171) . and 
the shapes are plotted in Fig. Eland HI The factorizable ansatz that is used to represent 
them in data analyses is given in Eq. (I6.22p and plotted in Fig. 

• Sinusoidal running and sharp feature. 

A sharp feature, in a potential or internal field space, introduces a sharp change in 
the slow-roll parameters, or the generalized slow-variation parameters. This can boost 
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the magnitudes of time-derivatives of some parameters by several orders of magnitude 
while still keep the power spectrum viable. These time-derivatives act as couplings in 
the interaction terms. So they enhance the non-Gaussianities among the modes which 
are near the horizon-exit. How deep they affect the modes inside the horizon depends 
on how sharp the changes are. 

The changes in these parameters can be roughly approximated as delta-functions in 
time. Correlation functions involve integrations of products of the slow- variation pa- 
rameters and the mode functions. The latter contain oscillatory behavior ~ e~ tKr , 
where the comoving momentum K is k\ + k 2 + k 3 for bispectra and r is the conformal 
time. The delta-function specifies a scale r*. This is why after integration the bispec- 
trum contains a sinusoidal factor ~ sin^/fc*), where = — 1/r* is the momentum 
of the mode that is near the horizon-exit at the time of the feature. So the most 
important property of this type of non-Gaussianity is this characteristic running. A 
numerical result of the running behavior is plotted in Fig. [7J An ansatz is given in 
Eq. (16T32|) and (IBTMjl . 

• Resonant running and periodic features. 

The periodic features do not have to be sharp. They introduce a small background 
oscillatory component in the slow-variation parameters. On the other hand, the mode 
functions are also oscillatory before they exit the horizon. Their frequencies are high 
when they lie deep inside the horizon and become lower as their wavelengths get 
stretched by the inflation. They are frozen after the wavelengths become comparable 
with the horizon size H~ x . This means that their frequencies continuously scan through 
the range from Mp (or some other large fundamental scale) to H. Therefore as long 
as the background oscillatory frequency uj satisfies H < uj < Mp, at some point 
during the evolution the small oscillatory component in the slow-variation parameters 
will resonant with the oscillatory behavior of the mode functions, and cause a large 
constructive contribution to the integration. 

The periodicity of the features leads to a periodic-scale-invariance in density pertur- 
bations. Namely, they are scale invariant if we rescale all momenta by a discrete e-fold 
27mH/u, where n is an integer. This is why the most important feature of this non- 
Gaussianity is a running behavior ~ sin(C\nK + phase), where C = u/H. This leads 
to the ansatz (16 .43 \\ . The full expression is given in (16.451) and plotted in Fig. [HJ A 
numerical result is plotted in Fig. [9j 

• Folded shape and non- Bunch- Davies vacuum. 

The usual mode function of the Bunch-Davies vacuum has the positive energy mode 
~ e~ tkr . Now we consider a non-Bunch-Davies vacuum by adding a small component 
of negative energy mode ~ e tkT . The three-point function involves an integration of 
the product of three mode functions with momentum fci, k 2 and k 3 . So the leading 
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correction to the Bunch-Davies results is to replace one of the fcj's with —ki. The usual 
K = k\ + &2 + ^3 in e~ %Kr becomes K = k\ + k% — k% and its cyclic. This effect is most 
important if factors of K appear in the denominators after the r-integration. Hence the 
most important feature of this type of modification is to enhance the non-Gaussianity 
in the folded triangle limit. An example of these bispectra is given in Eq. (16. 50}) and 
plotted in Fig. [101 Ansatz that partially capture this feature are given in Eq. (16.511) 
and (I6.52p , and plotted in Fig. [TTJ 

• Intermediate shapes and massive isocurvatons. 

All mechanisms discussed so for single field inflation apply to multifield inflation. We 
now consider new effects caused by introducing more fields to inflation models. These 
extra fields are called isocurvatons. 

Since light fields typically acquire a mass of order H, the Hubble parameter, we first 
consider the quasi-single field inflation models where there is one massless inflaton 
while the isocurvatons have mass of order H instead of massless. 

Unlike multifield slow-roll inflation, where each flat direction only has small non-linear 
terms in order to satisfy the slow-roll conditions, massive directions are not inflation- 
ary direction and are free to have large non-linear self-interactions. These non-linear 
interactions can be transferred to the curvature mode through couplings and source 
the large non-Gaussianity. 

The massive isocurvaton eventually decays after horizon exit simply because they are 
diluted by the expansion. How fast it decays depends on its mass. If the mass is heavier, 
m > \/2H, it decays faster. So the interactions can only happen when all modes are 
all closer to the horizon exit. This is closer to the case of the equilateral shape that 
we encountered above, and results in bispectra with quasi-equilateral shapes. If the 
mass is lighter, m < y/2H, it decays slower. More non-Gaussianity is generated in the 
super-horizon scales. This is closer to the case of the local shape that we will come to 
below, and results in bispectra with quasi-local shapes. Overall, at the squeezed limit 
^3 "C k\ = &2, the bispectrum shapes behave as {k^/ki) l / 2 ~ u , where v goes from to 
3/2 (corresponding to m from 3H/2 to 0) in the example we studied. In particular, 
if we take the massless limit while keeping the cubic self-interactions of isocurvaton 
large, we get a large bispectrum that has the same squeezed limit shape as the local 
one. Therefore, we have a one-parameter family of shapes that lie between the local 
and equilateral shape. 

The numerical results of these shapes are presented in Fig. [TU A simple ansatz is given 
in Eq. (17.181) that represents this family of shapes quite well, and is plotted in Fig. [T51 

• Local shape and massless isocurvatons. 

The fluctuation amplitudes of massless scalars do not decay after the horizon exit, 
and therefore this opens up a multifield space for the superhorizon evolution. For 
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superhorizon modes, we can use the separate universe picture and study the classical 
behavior of different patches of universe. These patches are separated by horizons and 
evolve independently of each other. So the evolution is local in space. 

Non-Gaussianities are generated when this multifield evolution is nonlinear, and any 
nonlinearity arising in the separate universe picture should also be local in space. A 
locality in position space translates to a non-locality in momentum space. This is 
why the resulted local shape bispectrum peaks at the squeezed limit. The behavior is 
{ks/ki)^ 1 for k 3 <C k\ = ki. This bispectrum is given in Eq. (18. 101) and the shape is 
plotted in Fig. [16j 

In all cases, the power spectra are either approximately scale-invariant so indistinguish- 
able from the simplest slow-roll models, or modified with features that can be made small 
enough to satisfy the current observational constraints. 

Large bispectra generically implies large trispectra, i.e. the four-point correlation func- 
tions. But trispectra contain more information and can be large even if bispectra are small. 
Experimentally, trispectra are more difficult to detect, but contain much more shape con- 
figurations. Each category above should have interesting extensions to trispectra. See 
Ref. |185H191j for the equilateral case and |192H197j for the local case. 

It is certain that this list will grow in future works, providing more refined and diverse 
connections between theories and experiments. 

9.2 A consistency condition 

As we have seen, in single field inflation, the mode that has exited the horizon is frozen. 
This is characterized by a constant ( over a horizon size patch. The physical meaning of the 
constant ( is a small rescaling of the scale factor. This is the only effect that the superhorizon 
mode has on modes with much shorter wavelength. This fact is used by Maldacena to derive 
a consistency condition (17] for the three-point correlation functions in the squeezed limit 
for single field inflation. 

• Consistency condition. In the squeezed limit k 3 <C k\ = k 2 , k 3 is the superhorizon 
mode that exited the horizon and acts as a zero-mode modulation to the two remaining 
modes. The correlation (Ck 1 (k 2 Ck 3 ) is an average of the following quantity 

(OiO^Cfa (9-1) 

over different (k 3 - We will shift the average (Cfc 3 ) to zero by definition. Here the two-point 
average (CfciCfc 2 )Ck 3 * s evaluated with different local scalings determined by the shift (k 3 . If 
the two-point function is exactly scale-invariant, (CfciOb)^ ^ s J us ^ a constant. So the 3pt 
vanishes because {(k 3 ) = 0. The non-zero contribution comes from the breaking of the scale- 
invariance. To see this, we expand the two-point average in terms of a long wavelength mode 
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( k4 near the scale (Cfe 4 ) = 0, 



(CfciCfc 2 )cfc 4 - (Cfc 1 Cfc 2 )o + 



^(CfciCfc 2 ) 



A'4 



(9.2) 



Multiply this with and average over C^. The first term contributes zero since this is the 
scale- invariant component. The second term gives 



^(CfciCfc 2 ) 



d(k 4 



(9.3) 



To get non-zero average, k 3 + k 4 = is needed. Using the relation d(k 4 = —din k, we get 



d((k 1 Ck 2 ) 



dink 



(Cfc 3 Cfc 3 ) • 



(9.4) 



The higher order terms in (I9.2p give 



2 (din A;) 2 



id 3 (c fel a 2 ) 



6 (dlnife)* 



<a 3 c fc 3 4 >+- 



(9.5) 



where k 3 + 2k 4 = 0, k 3 + 3k 4 = 0, and so on have to be satisfied respectively for each term 
to get non-zero average. If we only consider the tree-level three-point function, these higher 
order terms can be truncated since they involve more factors of P^ and should be related to 
the loop diagram contributions to (CfciCfc 2 Cfc 3 )- The tree diagram is 0(P£). 

To connect the averages we used here with the correlation functions that we defined 
in previous sections, we need to restore the phase factors. Here the two-point average 
(CfciCfc 2 )here is performed with the special point k x = k 2 in the phase space. To connect 
this with the previous definition of (CfciCfe)) we nee d to include the phase space in the 
neighborhood. Namely, (CfeiCfea) = (CfciCfc 2 )hcrc(2vr) 3 5 3 (k 1 + k 2 ), so from the definition (12.26P 
we have (Cohere = (2tt) 2 P c (A; 1 )/2A; 3 . Similarly, (C 3 ) = (C 3 )here(27r) 3 5 3 (£ k,). With the 
usual definition of the spectrum index n s — 1 = dlnP^/dlnk, from (I9.4p we get the following 
consistency condition (47] 



(CmCM -K - 1 )^p|n(^i)n(^)(2vr) 7 5 3 (^k t 



(9.6) 



Although originally derived for slow-roll inflation, the only assumption is the single field. 
So this applies to any single field inflation models and has important physical implications 
that we discuss shortly [198] . Note that the derivation of this relation (19.61) does not rely on 
the smallness of the slow- variation parameters either. For the general single field inflation 
models that we studied in Sec. 16.11 at tree level this has been checked with explicit results to 
three different orders [199( 150] including the slow- roll limit |47] . For resonance models, this 
is checked to the leading order [111] . 
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There are three types of interesting corrections to the condition (19. 6p . 

Firstly, as mentioned, the right hand side of ( 19.61) should receive corrections from loop 
contributions. These loop contributions are associated with higher derivatives of the two- 
point function. The terms (I9.5p . together with (19 .4p . provide the corresponding consistency 
conditions at different orders of P^. Note that for such orders, the correlation functions such 
as (CfciCfc 2 ) on the right hand side should also include loop corrections. 

Secondly, when we assume that the only effect of the frozen superhorizon mode on the 
much shorter scale is a constant background rescaling, we assume that there is no interaction 
when these modes are all within the horizon^ However, large subhorizon interaction is 
possible in some cases, such as in Sec. 16.31 and 16.41 Such interactions disappear below a new 
length scale at subhorizon, only then the above assumption becomes valid. For example 
in the resonance model, for H/u < k 3 /ki <C 1, the modes k 1; k 2 and k 3 are guaranteed 
to resonant with the oscillatory background at some point when all of them are within the 
horizon. Only if k 3 /ki <C H/u, such resonance will not happen. This is why the consistency 
condition is satisfied only in the very squeezed limit (k 3 <C k\jC with C = u/H ^> 1). For 
the squeezed region k\jC < k 3 -C fcj, the left hand side of the condition is larger than the 
right hand side by a factor of Ck 3 /2k\. For the non-Bunch-Davies vacuum case, a similar 
scale is determined by the UV cutoff r c , from which the non-Bunch-Davies vacuum starts to 
take effect. 

Thirdly, even after the long wavelength mode exits the horizon, as long as k 3 jk\ is not 
infinitely small, there is still dependence of the two-point function on the derivative of this 
long wavelength mode, in addition to the overall constant shift. This introduces a different 
type of finite k 3 jk\ corrections. They start from the second order in k 3 /ki, because the first 
order corresponds to the first spatial derivative of the long wavelength mode and should 
vanish due to isotropy [198J . These corrections will be amplified by the associated amplitude 
/atl" 100 ' an d gi ye an additive correction f™l~ loc {k 3 /ki) 2 to the n s — 1 on the right hand side 
of the condition. For a large /^ _loc , therefore, the condition needs to be satisfied in a very 
squeezed limit. The equilateral bispectra (16 . 1 6[) and (I6.17P are this type of examples. 

The consistency condition ( 19.61) can be straightforwardly generalized to higher order 
correlation functions |192[ll87j . We emphasize that this condition only applies to single field 
inflation models. For inflation models involving more than one field, as we have seen, non- 
Gaussianities can be transferred from the isocurvature directions which do not respect this 
relation. 

• Physical implication. Besides providing consistency checks for analytical computa- 
tions, the condition also has interesting physical implications. In the following, we discuss 
the scale invariant cases |198j , as well as the feature cases and loop corrections, ending with 
some cautionary remarks. 

This consistency relation implies that the tree-level bispectrum in the squeezed limit is 
determined by the power spectrum and spectral index. We distinguish the following two 

2 I would like to thank Yi Wang for helpful discussions on this point. 
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cases. For the scale-invariant case, n s — 1 is of order 0(e) and the right hand side of (19. 6p 
takes the local form. Indeed, as we have seen, for the single field inflation models where the 
non-Gaussianities are large, they take the equilateral forms which vanish in the infinitely 
squeezed limit. For the non-scale- invariant case, especially the highly oscillatory case such 
as the resonance model, the power spectrum can be highly oscillatory and n s — 1 becomes 
large. This can still be consistent with observations since the large n s — 1 is also highly 
oscillatory and therefore may escape a detection so far. But such a running non-Gaussianity 
is orthogonal to the scale-invariant forms. 

For the loop diagrams, in the scale-invariant case, these terms are suppressed by higher 
orders of slow- variation parameters from, e.g. d 2 {C,k 1 Ck 2 ) / {d In k) 2 , and higher orders of ( from, 
e.g. (C| 3 ); in the non-scale-invariant example, the extra terms are still highly oscillatory. 

In summary, a detection of an approximately scale-invariant local non-Gaussianity in 
the infinitely squeezed triangle limit with > 0(e) can rule out all single field inflation 
models. 

In experiments, however, the triangle cannot be perfectly squeezed. So it is an important 
question how squeezed it should be to achieve the above goal. For example, in the third 
type of corrections we discussed previously in this subsection, we need /a^ - loc (^3/ ^i) 2 to 
be smaller than n s — 1 for the consistency condition to hold, so that the contaminations 
from whatever non-local f^L is small. Assuming the primordial local form is practically 
detectable only if f°l > 0(1), we at least need fNL~ loc (h/ki) 2 < 0(1). For the class of the 
general single field models we studied in Sec. El if the other forms of non-Gaussianities, such 
as the equilateral one, can be constrained below f™i~ loc ~ 0(10), a squeezed configuration 
with k^/ki < 0.1 will be enough for our purpose. However, a completely model-independent 
statement is much trickier, because there may be bispectra with very large amplitude but 
orthogonal to any known bispectra that have been constrained experimentally. Besides 
that, in the second type of corrections, large finite-fcs/fci corrections can also arise due to 
subhorizon interactions. Therefore, as a cautionary remark, if we would like to rule out 
all single field inflation models in a rigorous model-independent fashion with a detection of 
scale-invariant local non-Gaussianity, we have to keep in mind the caveat that there may be 
single field models which only respect the consistency condition in a very squeezed region 
beyond the experimental reach. 

9.3 Superpositions 

Different shapes and runnings of non-Gaussianities can be superimposed in inflation models. 
For example, 

• Mixing shapes. 

It is possible that different non-Gaussianity generation mechanisms are from different 
components in a model, or at different stages during inflation. So two or more different 
shapes can get mixed, and the final shape can be rather different. For example, in 
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1.0 00 

Figure 17: A mixing of the equilateral (Fig. H]) and local shape (Fig. IT6i) . 

Fig. [13 we plot a mixed shape between the local and equilateral shape. Notice that this 
is different from the intermediate shapes, since obviously the squeezed limit is always 
dominated by the local form. Examples of such models are discussed in Ref. |200] . 

• Mixing shape and running. 

The shapes can also be mixed with runnings. Same as the power spectrum, the non- 
Gaussianities generically have some mild scale dependence. But a more dramatic 
case is the superposition with a strong running, such as the sinusoidal or resonant 
running. For example, an inflaton passing through features frequently and turning 
constantly at the same time on a potential landscape can generate a bispectrum which 
is a superposition of the resonant running and intermediate shape, as we illustrate in 
Fig. HHJ Clearly, these two signals are orthogonal to each other very well, and have to 
be picked up separately through different methods in data analyses. 

• Orthogonalization. 

If a non-Gaussianity is the linear superposition of several base components, one can 
generally perform a change of bases to make the new bases orthogonalized. For ex- 
ample, as we have seen in Sec. I6.1[ the leading large bispectrum has two components, 
S\ and S c . The two shapes are very similar, and represented by the equilateral ansatz 
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Figure 18: A mixing of an intermediate shape [u = 7/6 in ( IT. 181) ] and a resonant running 

(ES]). 



in data analyses. However since they do have small difference, one can subtract their 
similarities and get a new orthogonalized base component |173j . The orthogonalization 
is defined by the shape correlator such as f)5.12p . Using this definition, the new bases 
can be chosen as 

Si « S A + 0.22S C (9.7) 

and S2 = S c . [Note that the S\ and S c used here do not include the pref actors 
(1/Cg — 1 — 2A/S) and {l/c 2 s — 1) in (I6.16P and (I6.17p .] Their shapes are shown in 
Fig. [19j Notice that Si is half positive and half negative. Because Si is not of the 
simplest factorizable type, the following simple ansatz has been proposed to represent 
Si in data analyses |173j . 

Ansatz = "IS (J^ + 2 perm.) + 18 f^- + 5 perm.) - 48 . (9.8) 

We plot the shape of this ansatz in Fig. [201 The current CMB constraint on this 
orthogonal ansatz is —410 < /$£ h < 6 []]. 

For known examples of general single field inflation, such as the DBI and k-inflation, 
we generically get equilateral shapes. This is also clear from their physical origin that 
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+ = + 




Figure 19: Orthogonalization of two shapes in Sec. 16.11 (Fig. |3]and Fig. HJ) through a change 
of bases, c\S\ + c c S c = C\S\ + 0262. 




Figure 20: An ansatz —S°^ tz in (19. 8 j) for the orthogonal shape Si in (19. 7ft . Note we added 
a minus sign in this plot. 
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Figure 21: Another factorizable orthogonal ansatz (19. lip . 

we have emphasized. The orthogonal shape relies on a delicate cancellation between 
the two generic shapes. In principle one can do this since the required parameter 
space is allowed in our effective field theory of general single field inflation in Sec. 16.11 
and this may provide guidance to future model building. For example, one may fine- 
tune the parameters in the k- inflation models [93, 94J . Therefore, unlike the previous 
cases, the direct motivation here is more oriented to data analyses. The advantage 
of this operation is that it makes full use of data, which impose constraints on both 
components. In addition, as a bonus, the ansatz for the equilateral ( 16. 22ft . folded 
(I6.5ip and orthogonal (19.81) shapes are not linearly independent. As we can see, they all 
happen to be the equilateral ansatz shifted by a constant shape ansatz (5* = const.) [63] . 
Constraining two orthogonal bases provide efficient constraints on all three of them. 

Let us do a more data-analysis oriented exercise. We would like to construct an ansatz 
that is orthogonal to both local and equilateral ansatz, since both were well constrained 
by data. (Note that S^^ tz in (I9.8P is not quite orthogonal to the local ansatz, with 
a correlation ~ —0.48). To do this we start with a trial shape St r iai> and demand the 
new orthogonal shape 

^ansatz = Atrial + CiS 1 " + C2»S' a n Satz (9-9) 

be orthogonal to both the local and equilateral ansatz, 

oorth,2 oloc _ r\ riorth,2 <-,eq n (Q ~\f\\ 

^ansatz J u i ^ansatz ' ^ansatz u ■ (V.l-UJ 

The simplest factorizable trial shapes can be either the constant shape or the local- 



72 



like shape ki/k 2 + 5 perm., and both give the same result. Let us use the constant 
shape Striai = 1 as an example. Solving the conditions (I9.10j) gives c\ = —0.0953 and 
C2 = —0.204. So the new orthogonal ansatz is 

S£H = 1.19 (JL + 2 perm.) - 1.22 f| + 5 perm.) + 3.44 . (9.11) 

The numerical details may change slightly depending on the detailed definition and 
computation of the inner product (15. lip . The shape is plotted in Fig. [23] It is somewhat 
exotic but the ansatz is simple. By construction, this ansatz is much more orthogonal 
to the local form (with correlation ~ 0) than the S°^ tz currently used in Ref. |173[ [Tj. 
It also happens to have reasonably large correlation (~ 0.86) with the orthogonal shape 
in single field inflation (the Si shown in Fig. [T^]) . similar to that (~ —0.91) between Si 
and Ansatz- Obviously, other choices of trial shapes can result in more exotic orthogonal 
shapes. 

One can perform a similar orthogonalization for the two shapes in ( I4.22p . now they 
are both local to start with. More generally, if a non-Gaussianity has more base 
components, we can orthonormalize all of them one by one, in the sense of the Gram- 
Schmidt process. 

9.4 Conclusion 

The field of primordial non-Gaussianity is growing rapidly in recent years, with simultaneous 
progress from the experimental results, data analyses methods, non-linear cosmology theo- 
ries, physical model buildings, computational techniques, and theoretical formalisms. The 
progress that we have seen so far is no doubt just a beginning. 

In this review, we have studied the primordial non-Gaussianities coming from the inflation 
models, especially various mechanisms that can produce observable large non-Gaussianities 
with viable power spectra. We emphasized the fingerprints that different underlying physics 
leave on non-Gaussian profiles, which break the degeneracy of model building. We described 
the physical pictures and presented their effective Lagrangians to the extent that they can be 
recognized when encountered in the inflation model building in a more fundamental theory. 
We also derived the resulting bispectra and represented them in terms of simple ansatz 
to the extent that they can be useful to data analyses. With the current rapid progress, 
we anticipate much more future developments along these lines through refinements and 
discoveries in both theories and experiments. 

The standard model of cosmology - the Big Bang theory with ACDM - is now estab- 
lished better than ever, with the precision data coming from the cosmic microwave back- 
ground and large scale structures. New data will continue to flow from many ongoing and 
forthcoming experiments. Although Nature does not seem to be obligated to provide us any 
more information beyond the standard model, exciting possibilities exist that would help 
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us to understand the origin of the Big Bang. These include the more detailed deviations 
from the scale-invariance of the power spectrum, the primordial gravitational waves that we 
may detect from the CMB polarization, the isocurvature perturbations between matter and 
radiation, and the primordial non-Gaussianities. Without these types of data, the number 
of theoretical models with degenerate observational consequences proliferate with time and 
it will be hard to understand the microscopic nature of the inflation beyond our current 
knowledge, as well as to distinguish inflation from other possible alternatives. As we have 
reviewed, primordial non-Gaussianity - the collider in the very early universe - is one of the 
few hopes. It is becoming a target of many modern experiments. We do not know which 
cards Nature is hiding from us, but we are hoping and preparing for the best. 
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